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PREFACE 

This  report  was  prepared  by  Darryl ).  Calkins,  David  S.  Deck,  and  Dr.  Devinder  S.  Sodhi,  Re¬ 
search  Hydraulic  Engineers,  of  the  Ice  Engineering  Research  Branch,  Experimental  Engineering 
Division,  U.S.  Army  Cold  Regions  Research  and  Engineering  Laboratory.  Funding  for  this  research 
was  provided  by  the  U.S.  Army  Engineer  District,  Detroit 

Dr.  Bruce  Pratte  of  the  National  Research  Council  of  Canada  and  Dr.  George  Ashton  of  CRREL 
technically  reviewed  the  manuscript  of  this  report 

This  project  has  involved  many  people  in  addition  to  the  authors  of  this  report  Calvin  Acker¬ 
man’s  skillful  and  persistent  craftsmanship  in  the  model  construction  made  the  work  go  very  smooth¬ 
ly.  He  also  participated  in  the  model  calibration  and  testing  programs.  The  tests  required  many 
people,  including  Christopher  Gernhard,  Carl  Martinson,  and  many  others  in  the  Ice  Engineering 
Research  Branch,  for  weighing,  measuring,  and  breaking  ice.  Winter  data  collection  was  also  made 
easier  with  helicopter  support  from  the  U.S.  Coast  Guard  Air  Station  in  Detroit  and  field  support 
from  the  Hydraulics  and  Hydrology  Section  of  the  Detroit  District. 

The  technical  assistance  received  from  Ronald  Wooley  and  other  personnel  of  the  Hydraulics 
Division  of  the  U.S.  Army  Engineer  Waterways  Experiment  Station  is  greatly  appreciated.  Their 
experience  and  expertise  in  model  calibration  made  the  job  go  very  smoothly.  A  review  panel, 
consisting  of  Samuel  Lazier,  Dr.  jean-Claude  Tatinclaux  and  Louis  Shows,  provided  the  authors  tech¬ 
nical  feedback  on  procedures  and  test  plans. 
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HYDRAULIC  MODEL  STUDY  OF  PORT  HURON 
ICE  CONTROL  STRUCTURE 


D.J.  Calkins,  D.S.  Deck  and  D.S.  Sodhi 


INTRODUCTION 

In  1978  the  U.S.  Army  Engineer  District,  Detroit,  requested  thatCRREL  undertake  a  study  to 
evaluate  various  methods  for  ice  control  at  the  outlet  of  Lake  Huron  because  of  potential  problems 
associated  with  extended  winter  navigation  in  this  area.  Several  studies  were  needed  to  properly 
determine  the  type  of  ice  control  structure  (ICS),  its  location  and  design  forces. 

This  report  focuses  on  the  results  of  a  refrigerated  hydraulic  model  study  for  determining  the 
optimum  siting  of  an  ICS.  Previous  modeling  studies  of  ice  control  structures  with  provisions  for 
navigational  openings  have  ali  been  conducted  in  warm  environments  using  synthetic  ice  (Calkins 
and  Ashton  1975,  Acres  American  1975,  Arctec  1978). 

The  Ice  Engineering  Facility  at  CRREL  allows  us  to  operate  a  physical  model  in  a  cold  environ¬ 
ment  where  natural  ice  can  be  grown  for  the  experimental  work.  In  addition,  synthetic  ice  can  be 
tested  for  comparison  with  the  data  from  the  natural  ice  tests.  In  this  study,  field  data  provided  by 
the  Detroit  District  from  the  performance  of  the  St.  Marys  River  ice  boom  (Acres  American  1975) 
could  be  compared  with  the  modeling  data. 


SCOPE  OF  WORK 

The  basic  objective  of  the  Lake  Huron  ice  control  study  was  to  provide  the  necessary  criteria  for 
detailed  design  of  the  ice  control  structure  at  Port  Huron.  The  two  major  design  criteria  for  the 
ICS  are  1)  the  configuration  and  2)  the  forces.  The  configuration  of  the  ICS  will  be  determined  by 
three  factors:  1)  orientation,  2)  location  and  3)  size  of  opening.  The  forces  on  the  ICS  are  com¬ 
posed  of  four  primary  components  acting  on  the  ice:  1)  wind,  2)  water,  3)  ship-induced  forces, 
and  4)  wave  forces. 

The  diversity  of  factors  affecting  the  two  design  criteria  was  such  that  the  research  was  divided 
into  9  studies,  some  requiring  modeling  and  others  not. 

The  following  factors  were  modeled  in  the  laboratory  to  determine  their  impact  on  the  ICS: 

1)  ice  discharge  through  different  size  openings,  2)  ice  forces  and  3)  orientation  locations. 

Six  other  studies  did  not  require  modeling:  1)  estimation  of  ice  discharge  for  natural  conditions, 

2)  analysis  of  wind  data  at  Port  Huron,  3)  evaluation  of  present  and  future  ice  control  structures 


subjected  to  ship  passage,  4)  impact  of  extending  the  ice  cover  on  Lake  Huron,  5)  physical  and 
photographic  documentation  of  ice  conditions  at  the  site,  and  6)  winter  measurement  of  flow  veloo 
ity  beneath  the  ice  cover. 


ICE  DISCHARGE  FROM  LAKE  HURON  INTO  ST.  CLAIR  RIVER 

Figure  1  is  a  general  location  map  of  the  major  water  and  land  masses  in  the  area  of  Port  Huron, 
Michigan.  The  average  annual  surface  ice  discharge  ( Decern ber-April)  passing  from  Lake  Huron 
into  the  St.  Clair  River  is  approximately  907  km2  (350  mile2)  as  gathered  from  five  years  of  record. 
Since  observations  were  made  only  during  daylight  hours,  this  value  is  twice  the  total  of  the  day¬ 
light  observations.  The  maximum  and  minimum  annual  values  are  1814  km2  (700  mile2)  and  438 
km2  (169  mile2).  The  average  value  of  the  ice  discharge  represents  less  than  1 .5%  of  the  surface 
area  of  Lake  Huron. 

The  ice  discharge  in  terms  of  the  percentage  of  surface  area  passing  through  the  305-m  (1000-ft, 
entrance  width  to  St.  Clair  River  was  abstracted  from  daylight  time-lapse  photography  for  six  wint< 
seasons,  1971-1976.  The  flow  velocity  through  this  305-m  section  varied  due  to  fluctuations  in 
lake  level  and  flow  discharge,  but  an  average  of  1.45  m/s  (4.75  ft/s)  was  used  to  determine  the 
surface  ice  discharge.  The  ice  thickness  could  not  be  measured  from  the  photography;  only  a 
relative  indication  of  thin  or  thick  ice  could  be  seen  by  observing  the  surface  conditions  of  the  ice 
floes  on  the  photos.  Figure  2  shows  airphotos  taken  by  the  Detroit  District  at  Port  Huron  of  the 
ice  discharge  on  23  January  1973. 

The  data  for  1975-76  contain  several  large  time  gaps  due  to  camera  malfunctions;  the  data  for 
the  other  five  years  have  only  a  few  missing  days  and  are  of  good  quality.  To  distinguish  between 
no  data  and  zero  ice  discharge,  a  value  of  2%  concentration  is  plotted  in  the  accompanying  figures 
to  represent  zero  ice  discharge. 

The  ice  concentrations  were  measured  by  projecting  the  time-lapse  photography  on  a  gridded  screen 
and  visually  averaging  the  daily  ice  runs  into  six  categories;  0,  0-10%,  11-25%,  26-50%,  51-75% 
and  76-100%. 


Figure  1.  General  location  of  study  area. 
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a.  Ice  discharge  ffom  Lake  Huron  into  the  St.  Clair  River. 


b.  Same  view  as  a  except  5  min  later. 


Figure  2.  Aerial  views  of  ice  discharge— Port  Huron,  23  Jan¬ 
uary  1973. 


The  ice  discharge  into  the  St  Clair  River,  in  percentage  concentration  for  the  months  of  March 
1972  and  March  1974,  is  shown  in  Figures  3  and  4  to  illustrate  the  variability  from  one  winter 
season  to  another.  The  within-season  variability  of  ice  discharge  can  also  be  seen  in  the  data  for 
January  and  February  1975  in  Figures  5  and  6.  Due  to  this  high  variability  from  day  to  day  and 
month  to  month,  a  seasonal  distribution  of  ice  discharge  is  needed,  without  regard  to  sequential 
events. 

Figure  7  is  a  graphical  representation  for  the  percentage  of  time  during  the  winter  season  that 
the  ice  discharge  is  greater  than  a  certain  value.  Three  values  of  ice  discharge  were  chosen:  13,19 
and  26  km2/day  (5.0,  7.3  and  10.0  mile2/day).  These  values  represent  surface  ice  concentrations 
at  the  entrance  to  the  river  of  35,  50  and  70%.  With  the  exception  of  the  1971-72  winter  season, 
the  number  of  times  large  quantities  of  ice  pass  from  Lake  Huron  into  the  St,  Clair  River  is  small, 
i.e.  less  than  1  day  in  10  and  more  likely  1  day  in  20. 

Approximately  15%  of  the  time  during  the  winter  season  the  ice  discharge  can  exceed  13  km2/ 
day,  and  only  5%  of  the  time  will  the  ice  discharge  exceed  26  km2/day,  excluding  the  abnormal 
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Fab  72 


Feb  '74 


Figure  3.  February  1972  Ice  discharge— con¬ 
centration  %. 


Figure  5.  January  1975  ice  discharge— concen¬ 
tration  %. 


Figure  4.  February  1974  ice  discharge-con¬ 
centration  %. 


Figure  6.  February  1975  ice  discharge— con¬ 
centration  %. 


season  of  1971-72.  During  1971-72  nearly  a  third  of  the  surface  ice  for  the  entire  year  flowed 
downriver  during  a  9-day  period  in  April.  All  other  years  had  very  minimal  or  no  ice  discharge  for 
April. 

A  correlation  of  wind  data  and  ice  discharge  at  Port  Huron  did  not  offer  any  predictive  relation¬ 
ship.  An  analysis  of  high  ice  discharge  rates  (11.5  km2  or  4.4  mile2/day)  for  a  minimum  of  two 
consecutive  days  showed  an  average  discharge  value  of  65  km2  (25  mile2)  per  event  before  the  ice 
re-arched  or  the  ice  concentration  became  too  low  for  the  ice  to  extend  upstream  of  Port  Huron  (25 
events  analyzed).  However,  the  standard  deviation  was  nearly  equal  to  the  mean  ice  discharge. 
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Figure  7.  Percentage  of  time  that  ice  discharge 
is  exceeded  on  any  day— Port  Huron. 


WATER  VELOCITY  PROFILES  AT  PORT  HURON 

The  velocity  distribution  of  water  flow  beneath  an  ice  cover  has  been  shown  to  closely  follow 
a  logarithmic  distribution  from  the  position  of  maximum  velocity  to  the  ice  cover  underside.  The 
gradient  of  the  velocity  with  depth  is  related  to  the  shear  stress  t,  acting  on  the  underside  of  the 
ice  cover.  The  velocity  distribution  function  is 

V\  =  2.5  V*  In  [^]  (1) 


where 

l/j  =  velocity  at  any  vertical  position  y 
V*  =  shear  or  friction  velocity  =  V Tj/pw 
pw  =  density  of  water 

ky  =  equivalent  sand  roughness  height  of  the  ice  cover  underside. 


Equation  2  can  be  integrated  over  the  region  affected  by  the  ice  cover.  The  lower  limit  is  fcj/30, 
where  the  velocity  was  assumed  to  be  zero  and  the  upper  limit  is  defined  as  Yt,  the  distance  to  the 
position  of  maximum  velocity  measured  from  the  bottom  of  the  ice  cover  (see  Fig.  8). 

The  mean  velocity  in  the  region  affected  by  the  ice  cover  is 


^mi  = 


2.5  V* 

”>T“ 


(2) 


which  can  be  further  simplified  by  neglecting  the  last  term  since  it  is  significantly  smaller  than 
the  other  two  terms  in  the  bracket: 


^mi  =  2.5  V,  (inj^Vl). 


(3) 


5 


Stationary  lea  Cover 


Figure  8.  Velocity  distribution  beneath  ice  cover. 


The  equation  for  the  maximum  velocity  is  thus 


''max  =  2-5  V,  [in )  ~lj  • 

Substituting  eq  4  into  3  yields 

i/  -  ^max  ~ 

2.5 

and  further  reduction  by  the  definition  of  the  shear  stress  as  pw  V2  yields 


/''max  ~  ''mi  \  2 

T|  Pw  \  ZS  )  ' 
Likewise  the  Darcy  friction  factor  f  is 


(4) 

(5) 

(6) 

(7) 


The  velocity  distribution  beneath  the  ice  cover  at  Port  Huron  was  measured  at  10  locations 
corresponding  to  previously  monitored  sites  during  open  water  conditions.  Table  1  gives  a  summary 
of  the  velocity  distribution  data  below  the  ice. 


Table  1.  Summary  of  velocity  profiles  on  Lake  Huron,  7  February  1979. 


Site 

Mean 

vel, 

(mis) 

Max 

vel. 

(m/s) 

Depth  at 
max  vel. 
(m) 

Ratio 

Y|/D, 

Friction 

velocity 

(m/s) 

Stress 

(N/m2) 

Friction 
coeff,  f 

Ice 

thick 

(m) 

36 

0.35 

0.44 

4.65 

0.47 

0.036 

1.27 

0.084 

0.22 

51 

0.27 

0.42 

5.15 

0.63 

0.001 

3.69 

0.411 

0.24 

49 

0.28 

0.35 

5.58 

0.78 

0.028 

0.81 

0.083 

0.22 

63 

0.28 

0.36 

5.15 

0.78 

0.034 

1.16 

0.122 

0.24 

32 

0.48 

0.60 

4.55 

0.41 

0.048 

2.34 

0.081 

0.24 

55 

0.45 

0.58 

4.55 

0.42 

0.051 

2.54 

0.099 

0.24 

30 

0,37 

0.46 

430 

0.69 

0.037 

1.38 

0.082 

0.29 

57 

0.30 

0.40 

5.62 

0.88 

0.039 

1.53 

0.134 

0.17 

61 

0.37 

0.45 

4.12 

0.56 

0.031 

0.97 

0.056 

0.26 

67 

0.36 

0.40 

4.59 

0.75 

0.016 

0.24 

0.015 

0.20 

The  basic  data  from  Table  1  relevant  to  this  study  are  the  shear  stress  values  of  the  ice  under¬ 
side.  The  mean  value  is  1.59  N/m2  (0.034  lbf/ft2),  with  a  standard  deviation  of  1.0  N/m2  (0.021 
lbf/ft2).  By  eliminating  the  highest  and  lowest  value,  Tj  =  1.5  ±  0.63  N/m2.  Other  useful  coeffic¬ 
ients  were  also  calculated  such  as  V *  and 

The  collection  of  the  field  data  was  a  time-consuming  and  slightly  dangerous  effort.  On  only 
one  day  in  15  during  three  separate  missions  to  Port  Huron  were  we  successful  in  making  the  neces¬ 
sary  measurements.  This  was  accomplished  by  airlifting  the  people,  velocity  probe,  ice  auger,  safety 
gear  and  life  raft  onto  the  ice  with  helicopters  (see  Fig.  9).  The  helicopter  support  was  provided  by 
the  U.S.  Coast  Guard  Air  Station  at  Detroit  Triangulation  data  for  the  velocity  measurement  sites 
were  radioed  to  the  pilots  by  field  personnel  from  the  Detroit  District,  Corps  of  Engineers, 
Hydraulics  and  Hydrology  Section,  who  set  up  on  the  horizontal  control  points. 

A  two-directional  magnetic  flow  meter  was  used  to  gather  the  data.  The  probe  is  a  3.75-cm-diam. 
sphere  with  a  resolution  to  2  cm/s  and  a  measuring  volume  of  1.5  sphere  diameters.  Velocity  meas¬ 
urements  were  attempted  for  the  full  depth,  but  the  drag  on  the  extension  rods  was  so  great  in  some 
locations  that  the  vibrations  in  the  rods  were  excessive.  After  bending  a  few  rods,  only  velocity 
measurements  to  the  depth  of  maximum  velocity  were  made. 


ICE  CONDITIONS 


The  mean  air  temperature  data  for  the  Port  Huron  area  were  abstracted  from  the  Great  Lakes 
ice  cover  atlas  distributed  by  thp  National  Oceanographic  and  Atmospheric  Administration  (Table 
2).  The  site  of  the  air  temperature  measurements  was  Alpena,  Michigan. 

The  following  table  summarizes  these  monthly  values  in  terms  of  degree  days.  The  solid  ice 
sheet  thickness  for  a  single  layer  can  be  calculated  by  assuming  a  freezeover  date  of  1  January  by 
using 


t  =  ov'S(°  days) 


Figure  9.  Measurements  on  the  ice— Port  Huron. 


Table  2.  Average  air  temperature  and  estimated  ice  thickness  for 
Lake  Huron,  Alpena,  Michigm,  data. 


Freezing  degree  days  Cumulative  Ice  thickness 


1.40  1.3 

L53  1.75 

1.58  1.9 

1.53  1.75 


.*** 

Jan 

277 

500 

277 

500 

Feb 

223 

400 

500 

900 

*v.  • 

Mar 

83 

150 

583 

10S0 

m 

April 

-83 

-ISO 

500 

900 

The  coefficient  a  must  account  for  the  local  geographic  setting,  snow  cover,  wind,  thermal  and 
other  meteorological  conditions.  Using  a  value  of  approximately  0.0238  m  day-5*  °C~^  (0.7  in.  day-^ 
°F“^),  the  maximum  thickness  is  roughly  0.58  m  (1.9  ft).  The  ice  thickness  on  Lake  Huron  can 
be  highly  variable  due  to  rafting  and  ridging,  with  thicknesses  of  1  m  (3  ft)  or  more  common.  The 
ice  thickness  distribution  near  the  shipping  track  can  also  be  highly  variable  and  thicker  due  to  the 
constant  breaking  of  ice  and  deposition  at  the  sides  of  the  track.  Large  floes  broken  by  vessel 
transits  were  observed  to  be  deposited  in  layers  of  3  to  4  floes  thick  on  each  side  of  die  shipping 
track  near  Port  Huron. 

Airphotos  are  usually  very  descriptive  in  their  presentation  of  the  ice  conditions  at  a  particular 
location.  Figures  lOa-f  include  photography  taken  during  typical  ice  conditions  at  Port  Huron. 

Figures  lOa-c,  looking  uplake,  show  the  ice  arched  from  shore  to  shore  with  open  water  extending 
only  a  few  hundred  to  a  few  thousand  meters  above  the  Blue  Water  International  Bridge.  The  open 
water  areas  are  difficult,  if  not  impossible,  to  predict,  depending  on  ice  conditions,  wind,  and  ship 
traffic.  A  vessel  track  can  be  seen  in  the  upper  portion  of  the  solid  ice,  extending  into  the  lake. 

The  ice  conditions  at  Port  Huron  range  from  the  presence  of  only  small  ice  floes  at  freezeup,  to 
the  solid  cover  or  individual  large  floes  during  the  middle  of  the  winter,  to  the  total  breakup  of  ice 
in  the  spring.  The  ice  along  both  shorelines  can  be  very  thick  due  to  the  pileup  conditions  from 
high  winds  as  well  as  from  the  forces  in  the  sustained  ice  arching  that  develops.  The  thickness  can 
be  as  much  as  6  to  7  m  (20  ft)  depending  on  the  specific  location.  An  airborne  radar  system  was 
employed  to  measure  the  ice  thickness  distribution  in  the  lake,  but  the  data  have  not  been  fully 
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Figure  1 0a-f  (cont ’d).  Various  views  of  the  stable  ice 
arch  between  the  shore  of  Port  Huron,  Michigan,  and 
Sarnia,  Ontario,  taken  during  1971-1978. 


reduced  or  interpreted  at  this  time.  The  ground  truth  observations  taken  during  these  radar  flights 
show  ice  thicknesses  ranging  from  0.15  m  to  over  1  m.  The  large  thicknesses  were  the  result  of 
rafting  and  piling  of  ice  floes,  leading  to  ice  covers  of  one  to  four  floes  in  thickness. 

The  shore  ice  ridges  were  at  least  3  m  (10  ft)  high  in  some  locations  during  1979  and  residents 
in  the  area  talked  of  ridges  6  to  9  m  (20  to  30  ft)  high.  Apparently  more  historical  studies  or 
future  observations  are  needed  on  shore  ice  pileup  and  its  potential  effects  on  the  shoreline's  en¬ 
vironment. 


PHYSICAL  MODEL 
Basis  for  selection 

A  physical  model  of  the  site  was  needed  to  document  the  following  parameters:  1 )  ice  dis¬ 
charge  through  an  opening  in  the  ICS  during  ship  transits,  2)  water  drag  on  the  ice  due  to  the  ICS, 
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3)  ship-induced  forces  on  the  ICS  and  4)  velocity  patterns  if  permanent  structures  are  planned. 

The  scale  of  the  Froude  model  chosen  was  1 :85  for  both  horizontal  and  vertical  dimensions. 

The  model  scale  selection  was  primarily  based  oh  the  Reynolds  number  being  greater  than  1000 
for  ice-covered  conditions  at  a  water  temperature  of  0°C  in  and  near  the  navigational  channel. 

This  particular  scale  of  1 :85  also  allowed  for  sufficient  surface  area  upstream  of  any  proposed 
structure  for  vessel  transit  through  the  model  structure.  The  Reynolds  number  criterion  is  thus 

VDt 

N  =  — 1  =  2000  (6) 

r  v 

where  V  =  mean  velocity  over  the  depth  D] 

Dj  =  depth  of  water  from  ice  underside  to  the  bed 
v  =  kinematic  viscosity  of  water. 


The  Reynolds  number  of  2000  accounts  for  the  floating  ice  cover  changing  the  hydraulic  radius 
or  hydraulic  mean  depth  to  approximately  £),/2.  The  kinematic  viscosity  of  water  near  0°C  is 
nearly  1.83  times  the  value  for  water  at  20°C. 

To  operate  the  hydraulic  model  at  T  =  20°C,  the  open  water  Reynolds  number  must  be  at  least 
2500  to  3000.  This  will  ensure  that  the  model  can  also  be  operated  at  T  =  0°C  with  an  ice  cover. 

A  Froude  model  governs  when  the  predominant  driving  force  is  gravity,  open  channel  flow. 
Principles  of  open  channel  flow  can  be  applied  to  a  floating  ice  sheet  by  adjusting  the  cross-sectional 
area,  hydraulic  radius  and  roughness  coefficient  Since  the  Froude  model  has  ensured  dynamic 
similarity,  kinematic  and  geometric  similarity  are  implied.  The  relationship  is  simply 


'W'Vfp 


VJV p 


|<Ufp)<^M,/2 


=  1 


(9) 


where 


Nf  =  Froude  number 
V  =  mean  velocity  vector 
h  =  characteristic  length,  taken  as  the  depth 
g  =  acceleration  of  gravity 
m,p  =  model  and  prototype  conditions  respectively. 


Without  geometric  distortion  there  are  no  differing  scale  factors  for  time,  acceleration,  etc.,  within 
the  vertical  and  horizontal  directions.  The  scaling  laws  for  the  various  items  such  as  the  geometries, 
mass  flow,  pressure,  etc.,  can  all  be  calculated  very  easily. 


Table  3.  Scaling  conversions. 

Ratio  Factor 


Length 

85 

*-r 

Area 

^2// 2 

7225 

<M2 

Volume 

L^il  3 
Lp/l-m 

6.14X105 

U-r>3 

Time 

V»m 

9.219 

(M1/2 

Density 

pplpm 

1 

i 

Water  velocity 

V'm 

9.219 

M/2 

Water  discharge 

<V<?m 

6.66X1 04 

ar>5/2 

Fluid  shear  stress 

VTm 

85 

Lr 

Forces 

FplF  m 

6.14X105 

<M3 

The  more  commonly  used  scaling  factors 
in  this  report  are  given  in  Table  3  where  the 
subscripts  m  and  p  are  previously  defined 
and  r  refers  to  their  ratio. 

Description 

The  undistorted  scale  of  1 :85  for  the 
fixed-bed  Froude  model  provided  accurate 
reproduction  of  the  hydraulic  conditions, 
flow  velocity  and  direction.  The  model  sim¬ 
ulated  an  area  extending  from  Lake  Huron’s 
entrance  to  the  St.  Clair  River  (Fig.  11), 
approximately  1.6  km  (1  mi)  uptake. 

The  hydrographic  survey  was  supplied 
by  the  Detroit  District.  Template  cross 
sections  normal  to  the  flow  were  established 
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o  Point  Velocity  Location* 
•  Gage  Location* 


Yacht  Club 


Ontario 

Canada 


Figure  11.  Model  limits,  measuring  positions 
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a.  Initial  layout  of  mode!  templates. 


b.  Partial  completion  of  sand  base  for  the  model. 


c.  Nearly  completed  template  and  base  foundation  of  sand. 


Figure  12.  Typical  scenes  during  mode!  construction. 
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d.  Placement  of  5  cm  (2  in.)  of  mortar  to  complete  the 
model  bed. 

Figure  1 2  (cont’d).  Typical  scenes  during  model  con¬ 
struction. 


i 


a.  Steel  reinforcement  mesh  placed 
into  the  mortar  shell. 


b.  A  Iternate  drying  and  curing  of 
mortar. 


Figure  13.  Final  phases  of  model  construction. 
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c.  View  of  personnel  coating  the  model  with  fiberglass 
resin. 


Figure  13  (cant'd). 


every  103  m  (340  ft)  to  provide  accurate  bed  reproduction.  Figures  12a-d  show  typical  scenes 
during  the  model  construction  phase.  In  addition  to  the  usual  sand  fill-mortar  skin  construction 
for  the  model,  other  measures  were  taken  to  protect  the  model  from  alternative  freeze-thaw  cycles 
required  in  the  testing  program.  The  mortar  was  reinforced  with  a  10-gauge  wire  mesh  and  the 
entire  surface  was  sealed  with  a  resin  over  a  fiberglass  mat  (Fig.  13a-c).  These  additional  methods 
prevented  significant  thermal  cracks  and  minimized  the  migration  of  moisture  through  the  mortar 
skin  and  subbase. 

Pumps  supplying  0  to  0.03  m3/s  (500 gpm)  and  0  to  0.06  m3/s  (1000  gpm)  fed  the  33-m- 
(1 10-ft-)  long  inlet  trough  equipped  with  a  partitioned  weir  for  flow  control.  Each  pump  line  con¬ 
tained  a  magnetic  flow  meter  for  discharge  monitoring.  The  downstream  flow  control  consisted 
of  a  3.58-m-  (1 1 ,75-ft-)  long  flapper  gate  located  at  the  end  of  the  model  area. 

The  water  levels  were  measured  at  four  locations  in  the  model,  with  the  Ft.  Gratiot  gauge  (Fig. 

1 1)  used  as  the  control  gauge  for  the  study.  The  other  water  level  measurements  were  conducted 
in  the  center  of  model  along  an  approximate  stream  line  which  coincided  reasonably  well  with  the 
navigation  channel  (see  Fig.  11).  These  midchannel  water  level  measurements  were  performed 
from  an  overhead  carriage. 

Instrumentation 

The  major  instrumentation  procedures  consisted  of  calibrating  11  strain-gauged  aluminum  cir¬ 
cular  bars  in  the  two  horizontal  directions  and  writing  the  computer  programs  to  collect  the  data 
on  the  forces  transmitted  to  these  bars.  The  cross  sensitivity  between  gauges  for  each  bar  was 
calculated  and  the  stored  data  on  the  force  levels  reflected  this  correction  for  the  cross  sensitivity. 
The  sensitivity  of  the  bars  for  recording  small  deflections  was  limited  to  a  0.01 -N  (0.0022-lbf) 
change  in  force.  The  stability  of  the  electronic  data  collection  system  (a  Neff  300  controlled  by  a 
Hewlett-Packard  9845  computer)  was  excellent;  0.05  N  (0.01 1  Ibf)  was  the  maximum  drift  over 
a  three-  to  four-hour  period  for  any  one  bar  direction,  and  usually  the  drift  less  than  0.02  N. 

The  floating  ICS  was  attached  to  the  1 1  vertically  mounted  instrumented  bars,  which  measured 
the  load  in  the  two  horizontal  directions.  These  bars  were  attached  to  a  square  tubing  frame  that 
was  in  turn  supported  by  four  vertical  supports  (Fig.  14).  The  four  vertical  supports  were  located 
downstream  of  the  ICS  and  did  not  interfere  with  the  force  measurements  or  movement  of  ice 
through  the  gap  opening. 
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Figure  14.  Typical  view  looking  west  of  the  mode!  ship 
in  a  plastic,  randomly  shaped,  fragmented  ice  field. 
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At  the  beginning  of  the  study  the  type  of  ice  control  structure  that  would  be  acceptable  to  the 
surrounding  environment  was  not  known.  To  construct  a  model  ICS  structure  at  a  scale  of  1 :85 
and  maintain  the  exact  geometric  replications  was  not  practical.  Since  no  particular  ICS  was  being 
tested  a  compromise  on  the  model  structure  to  be  placed  was  made. 

A  series  of  floating  polyethylene  blocks  2.5  cm  x  7.5  cm  x  1.25  cm  thick  (1  x3x  0.5  in.)  tied 
together  with  nylon  line  were  held  in  place  by  vertical  rods.  This  floating  ICS  held  the  ice  sheet  in 
place  and  transmitted  the  loads  to  the  vertical  rods.  By  moving  a  frame  system  from  which  the 
vertical  rods  were  suspended,  the  ICS  could  be  moved  in  the  model.  The  ICS  was  divided  into  two 
sections,  each  beginning  at  the  U.S.  and  Canadian  shorelines  and  extending  toward  the  navigation 
channel.  Adjustments  in  the  support  frame  allowed  for  different  gap  openings  to  be  studied  in  the 
ICS. 

Mean  velocity  data  of  the  flow  beneath  the  ice  cover  at  selected  locations  were  collected  using  a 
0.95-cm-  (38-in.-)  diameter  two-directional  magnetic  flow  velocity  probe.  Reliable  documen¬ 
tation  of  the  vertical  velocity  profile  could  not  be  measured  in  the  model  with  this  probe  because 
it  was  too  large  in  comparison  to  the  water  depths  of  6  to  12  cm  (3  to  5  in.)  in  the  model.  There¬ 
fore,  velocity  measurements  were  made  only  to  estimate  the  mean  velocity.  Accuracy  is  ±  1.0 
cm/s  in  the  model. 

Open  water  calibration 

The  model  was  adjusted  to  simulate  conditions  determined  from  prototype  drogue  studies  that 
were  conducted  in  May  and  October  1978  by  the  Detroit  District  (Fig.  15).  Cylindrical  2-cm- 
(3/4-in.-)  diam.  PVC  floats  were  used  to  measure  the  streamlines  and  velocities.  The  floats  were 
adjusted  to  the  proper  corresponding  prototype  depths  of  measurement. 

The  various  partitions  of  the  inlet  weir  were  adjusted  until  the  flow  lines  in  the  model  were  in 
satisfactory  agreement  with  the  prototype  velocities,  both  in  magnitude  and  direction.  The  open 
water  calibration  was  performed  in  conjunction  with  personnel  from  the  Hydraulics  Division  of  the 
U.S.  Army  Waterways  Experiment  Station,  Vicksburg,  Mississippi. 

The  flow  used  for  model  calibration  (Fig.  1 5)  was  a  result  of  the  two  drogue  studies.  As  both 
flows  were  within  2%  of  each  other,  a  flow  of  5700  m3/s  (201,000  cfs)  with  a  Ft.  Gratiot  level  of 
578.87  ft  was  used.  Although  the  flow  lines  in  the  model  did  not  change  appreciably  for  different 
model  flows,  a  set  of  drogue  data  at  a  lower  flow  should  have  been  gathered  to  verify  the  model 
results.  Flows  at  different  gauge  elevations  did  not  significantly  change  their  velocity  patterns  or 
magnitudes  by  more  than  10%  over  the  entire  model. 


16 


Open  water  tests 

The  mean  annual  winter  flow  (December-April)  into  the  St.  Clair  River  was  calculated  to  be 
4704  m3/s  (166,000  ft3 /s)  with  a  standard  deviation  of  708  m3/s  (25,000  ft3/s).  After  consulta¬ 
tion  with  the  Detroit  District,  three  flows  were  evaluated:  4700  m3/s  (166,000  cfs),  5385  m3/s 
(190,000  cfs),  and  5950  m3/s  (210,000  cfs)  in  the  model.  The  largest  flow  represents  approximately 
the  maximum  monthly  flow  recorded  during  the  last  20  years  for  the  winter  months.  The  corre¬ 
sponding  water  surface  elevations  at  the  new  Fort  Gratiot  gauge  for  these  three  flows  are  1 75.62  m 
(576.19  ft),  176.22  m  (578.16  ft)  and  176.73  (579.8  ft),  International  Great  Lakes  Datum  ( IG ID). 
The  water  level  elevations  for  the  same  flow  will  vary  considerably  from  year  to  year,  and  a  0.7-m 
(2.3-ft)  difference  in  height  was  noticed  for  the  mean  winter  flow. 

The  mean  velocities  from  the  model  drogue  study  over  the  entire  model  for  the  three  flows 
mentioned  earlier  are  given  in  Figures  16-18.  The  general  streamline  patterns  are  the  same  and  the 
velocities  do  not  differ  significantly  in  the  center  of  the  model,  i.e.  the  navigational  channel.  A 
significant  change  in  the  flow  velocities  occurs  along  the  U.S.  shore  due  to  the  shallow  depths  ex¬ 
tending  into  the  lake  about  100  to  200  m  (328  to  656  ft). 


Experimental  procedures  and  techniques 

Both  plastic  ice  floes  and  natural  ice  grown  under  the  refrigerated  conditions  in  the  testing  area 
were  used  in  the  experiment.  The  testing  techniques  for  the  laboratory-grown  ice  were  different 
from  those  used  for  the  plastic  ice. 

The  ice  discharge  was  measured  by  two  methods.  If  plastic  ice  were  being  used,  then  one  merely 
had  to  count  the  number  of  floes  passing  through  the  opening.  The  tests  with  the  plastic  ice  used 
either  square  floes  of  uniform  size  or  randomly  shaped  floes  whose  average  size  was  measured  pre¬ 
viously  from  the  sample  population.  In  this  manner  the  area  per  floe  was  known  and  multiplication 
by  the  number  of  floes  released  gave  the  area  of  ice  discharged.  The  thickness  of  the  plastic  "ice” 
was  0.635  cm  (1/4  in.)  or  0.54  m  (1.77  ft)  for  the  prototype. 

The  measurement  of  the  laboratory -grown  ice  discharge  was  accomplished  by  weighing  the 
amount  of  ice,  determining  the  average  floe  thickness  and  counting  the  number  of  floes.  From 
the  weight  measurement  and  ice  thickness,  the  surface  area  would  be  computed.  The  average  floe 
size  was  calculated  by  counting  the  number  of  floes  per  disturbance  and  dividing  the  total  surface 
area  by  the  number  of  floes.  The  data  from  several  disturbances  (200-500  total  floes)  would  be 
used  to  calculate  the  mean  floe  size.  These  relationships  are 


A  = 

r  Pi'i 


(10) 


a=  [/ A  tIN ,]*  (11) 

where  A  r  =  area  of  ice  released 
W-  =  mass  of  ice  released 
Pj  =  density  of  ice  =  0.92 
fj  =  thickness  of  ice 
N j  =  number  of  ice  blocks  released 
a  =  average  fioe  size. 

Flow  establishment 

Once  the  flow  distribution  was  set  with  the  partitioned  weirs,  no  changes  were  introduced  at 
the  upstream  end.  The  water  level  at  the  downstream  portion  of  the  model  was  then  controlled 
by  a  gate  to  achieve  the  appropriate  prototype  values  at  the  Fort  Gratiot  gauge. 
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Figure  16.  Mean  velocities  for  4700  m3 /s— drogue  method- 
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Figure  18.  Mean  velocities  for  5950  m3 /s-, drogue  method. 


The  water  surface  slope  from  points  B  to  C  (see  Fig.  11)  was  measured  for  many  of  the  tests 
using  both  plastic  and  natural  ice,  and  little  difference  in  slope  was  detected.  In  fact,  measurement 
error  was  more  of  «.  roblem,  as  the  slopes  were  very  mild  (5  x  10-5)  with  a  vertical  drop  of  only 
0.4  mm  over  a  distance  of  8.534  m  (28.0  ft).  The  accuracy  of  the  point  gauges  was  only  0.1  mm. 

Artificial  ice 

The  placement  of  the  plastic  "ice”  on  the  water  surface  was  done  manually.  The  water  flow  in 
the  model  was  reduced,  and  then  the  flow  was  gradually  increased  to  the  desired  conditions.  All 
experiments  were  started  with  the  single  layer  ice  floes.  A  small  wooden  stick  was  placed  across  the 
opening  between  the  two  ice  control  structures  to  prohibit  any  ice  discharge  until  the  flow  condi¬ 
tions  and  water  elevations  were  established  correctly. 

Natural  Ice 

The  procedure  used  to  break  the  freshwater  ice  was  again  manual.  The  water  flow  was  either 
reduced  to  0.006  m3/s  (100  gpm)  or  completely  stopped,  and  the  desired  ice  thickness  was  achieved 
by  cooling  the  room  overnight  at  the  appropriate  temperatures.  Maintaining  a  small  flow  minimized 
the  ice  buildup  in  the  inlet  troughs  and  on  the  outlet  gate  and  produced  only  a  small  zone  in  the 
model  where  the  ice  thickness  was  not  uniform.  Otherwise  the  thickness  over  the  entire  ice  sheet 
was  constant  within  ±  1  mm. 

The  ice  was  first  broken  away  from  the  shorelines,  the  inlet  trough  and  the  gate  so  that  the  ice 
cover  would  float  when  the  water  discharge  and  levels  were  increased.  With  sufficient  flow  beneath 
the  ice  cover,  all  ice  downstream  of  the  ICS  was  broken  up  and  allowed  to  float  over  the  gate. 

During  this  time  the  room  was  being  warmed  to  an  air  temperature  slightly  above  freezing,  which 
was  necessary  to  prevent  refreezing  of  the  individual  ice  floes  broken  upstream  of  the  ICS. 

Since  the  study  did  not  concern  itself  with  modeling  the  breaking  of  a  solid  ice  sheet  by  a  ship, 
the  freshwater  ice  was  not  modeled  for  its  flexural  strength.  When  the  room  air  temperature  was 
above  freezing,  the  water  discharge  was  stopped  and  the  ice  floes  were  broken  into  randomly  shaped 
pieces  by  walking  carefully  over  the  ice,  and  further  breaking  them  into  the  desired  flow  size  with 
blunt  wooden  probes.  The  smallest  average  floe  size  for  any  test  was  10  cm  (model  size).  When  the 
ice  breaking  was  completed,  the  initial  readings  for  the  force  measurements  would  be  taken  and 
then  the  water  flow  would  be  increased  to  the  proper  conditions.  The  purpose  of  decreasing  the 
flow  was  to  minimize  the  layering  of  the  ice  floes  during  breakage.  As  with  the  plastic  floes,  a 
single  layer  of  fragmented  ice  was  the  starting  point  for  all  tests. 

The  model  ship  was  loaned  to  CRREL  by  the  Department  of  Marine  Architecture,  University 
of  Michigan.  Its  scaled-up  prototype  dimensions  are  310  m  (1 020  ft)  in  length  with  a  beam  of  29 
m  (95  ft)  and  a  draft  of  8.23  m  (27  ft).  The  ship  was  guided  by  a  pulley  system  and  driven  by  a 
variable  speed  motor.  To  eliminate  the  problem  of  turning  the  ship  around  after  every  transit,  the 
stern  section  was  replaced  by  an  identical  bow  shape.  This  allowed  us  to  run  alternate  vessel  pas¬ 
sages  up  and  down  through  the  ice  control  structure  and  ice  field  with  minimal  delay  and  without 
disturbing  the  ice  field. 

A  typical  test  procedure  would  generally  follow  the  guidelines  below: 

1.  Select  ICS  orientation,  gap  opening  and  establish  a  low  flow. 

2.  Grow  ice  sheet  or  place  plastic  “ice”  on  water  surface. 

3.  Eliminate  all  ice  downstream  of  ICS. 

4.  Break  ice  into  desired  flow  size. 

5.  Take  initial  force  reading  on  ICS  at  very  low  or  no  flow. 

6.  Adjust  flow  and  water  level  at  Ft.  Gratiot  to  desired  conditions. 

7.  Transit  vessel  uplake  through  ICS  opening. 

8.  Measure  ice  discharge. 

9.  Take  force  measurements  after  ice  has  stopped  moving. 

10.  Transit  vessel  downlake  through  ICS  opening. 

1 1 .  Measure  ice  discharge. 
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12.  Take  force  measurements  after  ice  has  stopped  moving. 

13.  Repeat  steps  7  through  12  until  ice  field  is  depleted. 

If  dynamic  force  measurements  were  being  taken,  a  continuous  measurement  of  the  load  on  all  1 1 
bars  was  monitored  during  the  ship  transit  through  the  gap  opening.  The  scan  rate  was  3  samples/ 
second  for  22  channels  of  force  data,  i.e.  1 1  bar  positions  measuring  forces  in  both  x  and  y  direc¬ 
tions. 

Ice  cover  calibration 

The  calibration  of  the  model  for  ice-covered  conditions  was  verified  by  setting  the  model’s  Ft. 
Gratiot  gauge  at  the  desired  winter  level  and  then  determining  if  a  properly  scaled  water  shear  stress 
was  achieved.  The  prototype  shear  stresses  calculated  from  velocity  profiles  beneath  the  ice  cover 
are  given  in  Table  1.  The  average  of  eight  measurements,  excluding  the  highest  and  lowest  readings, 
was  1.5  N/m2  (0.03  lbf/ft2),  s.d.  =  0.63  N/m2  (0.14  lbf/ft2)  with  a  range  of  0.81  to  2.54  N/m2. 

The  shear  stresses  of  the  water  against  the  ice  cover  in  the  model  were  calculated  directly  by 
isolating  the  ice  cover  from  the  shoreline,  determining  the  load  to  the  ice  control  structure  from  the 
instrumented  bars,  and  dividing  by  the  surface  area  of  the  ice.  Then,  the  stream-wise  component 
of  the  ice  weight  was  subtracted  for  all  tests  and  its  value  was  found  to  range  from  30  to  165  g 
(0.66  Ibf  to  0.363  Ibf)  depending  on  the  total  ice  mass  on  the  water  surface  slope  of  5  x  10"5. 

The  total  measured  load  varied  from  150  to  400  g  depending  upon  the  shear  stress  and  total  ice  mass. 

Two  techniques  of  isolating  the  ice  cover  were  evaluated  and  both  gave  good  results.  The  first 
technique  was  used  with  the  plastic  “ice.”  The  instrumented  and  floating  ICS  prevented  the  ice 
from  moving  downstream,  while  another  floating  ice  boom  surrounded  the  ice  floes  upstream  and 
was  connected  to  the  outside  ICS  force  measuring  rods.  This  outer  boom  prevented  the  ice  from 
coming  in  contact  with  the  shoreline.  The  shear  stress  measurements  for  a  flow  of  5950  m3/s 
(210,000  cfs)  were  reproduced  with  good  results  that  scaled  up  to  prototype  values.  The  first  test 
with  the  U.S.  portion  of  the  ICS  in  position  1  (see  Fig.  19)  had  four  separate  measurements  taken 
for  the  shear  stress  and  produced  a  scaled-up  value  of  1.28  N/m2  (0.027  lbf/ft2).  The  second  test, 
repeated  four  days  later  but  on  position  2  (Fig.  19)  gave  a  value  of  1.29  N/m2  (0.027  lbf/ft2).  The 
shear  stress  for  a  flow  of  4700  m3/s  (166,000  cfs)  at  578.1  ft  was  also  measured.  It  gave  a  value 
of  roughly  1.81  N/m2  (0.038  lbf/ft2). 

The  second  technique  used  a  solid  ice  sheet  covering  only  a  portion  of  the  model.  The  ice  cover 
shear  stress  for  a  solid  ice  sheet  could  be  measured  without  attaching  a  surrounding  boom  by 
physically  moving  the  large  sheet  (~  6  x  6  m)  to  different  sections  of  the  model  and  letting  the 
floe  come  to  rest  against  the  ICS.  The  shear  stress  computed  was  1.38  N/m2  (0.029  lbf/ft2). 

The  floe  extended  from  the  Canadian  ICS  to  range  4  in  the  model. 

Moving  this  floe  to  the  center  of  the  navigation  track  with  its  boundaries  between  ranges  22 
and  17  resulted  in  a  shear  stress  of  2.23  N/m2  (0.047  lbf/ft2),  a  value  very  consistent  with  the  equiv¬ 
alent  prototype  measured  conditions  at  sites  32  and  55.  The  values  of  the  shear  stress  are  high  be¬ 
cause  the  ice  floe  being  used  had  remained  after  a  series  of  tests  and  the  underside  had  roughened 
due  to  melting  caused  by  warm  water  at  0.4°C  entering  the  inlet  trough. 

Once  the  readings  were  completed  for  determining  the  shear  stress  values,  the  upstream  ice  boom 
surrounding  the  ice  floes  was  removed  and  the  ice  was  permitted  to  come  into  contact  with  the 
shoreline  so  that  force  readings  could  be  obtained  in  this  situation.  To  ensure  that  the  ice  had 
sufficient  contact  with  the  shoreline,  the  ice  cover  was  surged  by  raising  and  lowering  the  down¬ 
stream  control  gate,  producing  a  small  wave  upstream  that  moved  the  ice  cover  about  1  m  (3.28 
ft).  When  the  disturbance  decayed,  the  ice  floes  would  be  rearranged  and  good  contact  would 
occur  along  the  shoreline.  The  average  load  to  the  ice  control  structure  with  the  fragmented  ice 
in  contact  with  the  shoreline  was  approximately  80%  of  the  ice-free  load,  with  a  range  of  73  to 
83%  for  four  replicate  tests. 

Point  velocity  measurements  were  conducted  under  the  ice  at  the  same  locations  where  open 
water  and  ice-covered  velocity  profiles  were  gathered.  Velocity  measurements  beneath  fragmented 
plastic  and  natural  ice  covers  and  a  uniform  consolidated  (nonbroken)  ice  cover  were  conducted. 
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Table  4.  Mean  velocity  measurements  beneath  ice  cover  of 
fragmented  square  plastic  floes. 


Flow  rote  Flow  rote 


Site* 

5950  m3/s  (210,000  ft3 Is) 
Velocity 

4  700  m3/s  (166,000  ft3 Is) 
Velocity 

(m/s) 

(ftls) 

(m/s) 

(ftls) 

30 

0.56 

1.8 

0.53 

1.7 

57 

0.59 

1.9 

0.68 

2.2 

32 

0.65 

2.1 

0.59 

1.9 

55 

0.62 

2.0 

0.68 

2.2 

34 

0.53 

1.7 

- 

- 

36 

- 

- 

0.37 

1.2 

51 

0.34 

1.1 

0.31 

1.0 

49 

0.34 

1.1 

0.31 

1.0 

63 

- 

- 

0.40 

1.3 

61 

0.31 

1.0 

0.28 

0.9 

26 

0.59 

1.9 

“ 

* 

*  See  Figure  1 1  for  location  of  the  sites. 


Table  5.  Mean  velocity  measurements  beneath  a  consolidated  ice  cover 
(1  m,  [3.3  ft]  thick)  and  a  fragmented  freshwater  ice  cover  (0.76  m, 
[2.5  ft]  thick). 


Fragmented 

Consolidated  Ice  cover  ice  cover 


5/fe* 

5950 

m3/s 

210  000  5385 

ft3 fs  m3/s 

190,000 
ft3  Is 

4700 

m3/s 

166,000 
ft3  /s 

4700 

m3/s 

166,000 
ft3  /s 

(mis) 

(ft/s) 

(m/s) 

Velocities 
(ft/s)  (mis) 

(ft/s) 

(m/s) 

(ft/s) 

30 

0.47 

1.5 

0.40 

1.3 

0.28 

0.9 

0.56 

1.8 

57 

0.65 

2.1 

0.65 

2.1 

0.56 

1.8 

0.56 

1.8 

32 

0.68 

2.2 

0.84 

2.7 

0.84 

2.7 

0.93 

3.0 

55 

0.78 

2.5 

0.84 

2.7 

0.74 

2.4 

0.93 

3.0 

34 

0.62 

2.0 

0.74 

2.4 

0.65 

2.1 

0.78 

2.: 

36 

0.47 

1.5 

0.59 

1.9 

0.56 

1.8 

0.59 

1.9 

51 

0.47 

1.5 

0.56 

1.8 

0.47 

1.5 

0.40 

1.3 

41 

0.47 

1.5 

0.43 

1.4 

0.43 

1.4 

0.43 

1.4 

49 

0.34 

1.1 

0.34 

1.1 

0.34 

1.1 

0.43 

1.4 

63 

0.43 

1.4 

0.47 

1.5 

0.34 

1.1 

0.43 

1.4 

61 

0.40 

1.3 

0.47 

1.5  • 

0.47 

1.5 

0.37 

1.2 

♦See  Figure  1 1  for  location  of  the  sites. 


The  plastic  fragmented  cover  provided  the  more  uniform  condition  of  single  layered  ice  floes  over 
the  entire  model. 

Table  4  is  a  summary  of  the  mean  velocities  under  a  plastic  fragmented  ice  cover.  The  two 
flow  rates  in  the  table  show  similar  velocity  conditions  at  the  same  sites,  which  are  within  the 
range  of  the  velocity  measurement  equipment.  Table  5  gives  velocity  data  taken  under  a  labora¬ 
tory  grown  ice  sheet  for  several  flow  rates  using  the  same  consolidated  ice  sheet.  In  general  the 
mean  velocity  conditions  are  similar  at  each  site.  The  velocity  patterns  (streamlines)  with  a  solid 
ice  cover  were  determined  by  injecting  dye  beneath  the  ice  cover  and  observing  its  path.  Figure 
20  gives  the  drogue  lines  for  a  flow  of  5950  m3/s  (21 0,000  ft3/s).  The  patterns  are  generally  sim¬ 
ilar  to  the  open  channel  flow  near  the  center  of  the  lake  (navigation  channel),  but  because  of  the 
ice  blockage  along  the  shallow  U.S.  shoreline,  the  flow  tends  to  be  slightly  diverted  toward  the 
navigation  track  by  approximately  46  m  (150  ft)  between  ranges  1 7  and  21 . 


The  velocity  conditions  changed  more  near  the  U.S.  shoreline  than  in  any  other  location  due  to 
the  thick  ice  cover  that  moved  the  flow  easterly  into  the  lake.  This  is  seen  by  the  increase  in  the 
mean  velocity  at  sites  30,  32  and  55  with  hardly  any  change  in  mean  velocities  from  the  navigational 
channel  to  the  Canadian  shoreline.  Figure  21  represents  the  stream  lines  beneath  a  fragmented  ice 
cover.  Other  than  the  shift  in  the  flow  along  the  U.S.  shoreline  between  ranges  17  and  21,  there 
is  only  a  small  change  in  the  flow  patterns  with  respect  to  a  fragmented  or  a  consolidated  ice  cover. 


ICE  CONTROL  STRUCTURE  ORIENTATION 


The  orientation  of  the  ice  control  structure  was  evaluated  using  square  plastic  ice  floes.  During 
these  tests  no  discharge  of  the  ice  was  allowed.  The  hydraulic  conditions  for  all  tests  were  a  flow 
rate  of  5950  m3/s  (210,000  cfs)  at  an  elevation  of  176.72  m  (579.80  ft  IGLD)  at  the  new  Ft. 
Gratiot  gauge.  The  priorities  for  analyzing  the  ICS  locations  were  its  1 )  safety  with  respect  to  nav¬ 
igation,  2)  effectiveness,  and  3)  minimal  cost 

The  navigation  channel  from  Lake  Huron  makes  an  alignment  shift  of  27°,  approximately  610  m 
(2000  ft)  above  the  entrance  to  the  St  Clair  River.  The  ICS  must  be  located  upstream  or  down¬ 
stream  from  the  turning  area  to  minimize  or  eliminate  any  problems  associated  with  negotiating 
this  turning  maneuver  during  ice  conditions.  Since  the  downstream  position  is  not  satisfactory 
because  of  high  flow  velocities,  the  location  of  the  ice  control  structure  above  the  turning  point  is 
the  only  practical  site. 

The  performance  of  the  ICS  in  retaining  individual  floating  ice  floes  is  governed  by  the  flow 
velocity,  water  depth  and  ice  floe  thickness  one  desires  to  contain.  A  summary  of  laboratory 
studies  on  stability  of  individual  ice  floes  was  provided  by  Ashton  (1974),  who  recommended  the 
following  equation: 


where  Vc  =  critical  velocity  upstream  of  cover 
tj  =  ice  thickness 
H  =  flow  depth 
Pj  =  density  of  ice 


If  the  actual  flow  velocity  is  less  than  VQ  for  the  given  parameters  of  tx  and  H,  then  juxtaposition 
of  individual  floes  will  occur  without  rafting. 

It  has  been  observed  in  the  field  that  an  ice  cover  can  form  with  thin  3-cm  (0.1 -ft)  ice  against 
flow  velocities  of  up  to  0.6  and  even  0.7  m/s  (2-2.3  ft/s).  This  implies  that  the  cover  forms  a  min¬ 
imum  thickness  and  progresses  upstream,  and  the  thickness  is  governed  by  a  stability  criterion  such 
as  the  one  suggested  by  Ashton  (1974).  At  the  beginning  of  the  ice  season,  containment  of  thin 
ice  floes  behind  an  ICS  wili  always  be  a  problem,  but  once  the  ice  floes  begin  to  develop  a  thickness 
of  5  to  10  cm  (0.1 5  to  0.3  ft)  and  attain  a  floe  surface  dimension  of  a  few  meters,  ice  cover  thick¬ 
ening  and  uplake  progression  behind  the  boom  will  begin.  In  all  model  tests,  the  ice  control  struc¬ 
ture  was  located  in  areas  with  flow  velocities  of  less  than  0.6  m/s  (2.0  ft/s)  and  the  cover  thickness 
was  always  greater  than  0.46  m  (1.5  ft). 

The  ICS  from  the  eastern  side  of  the  navigation  channel  to  the  Canadian  shoreline  is  sited  above 
a  breakwater  used  to  protect  a  small  marina.  The  velocities  in  this  area  are  generally  low  [<  0.4 
m/s  (1.25  ft/s)  ] .  Also,  this  location  should  minimize  problems  with  near-shore  structures  and  not 
affect  the  marina  operation  during  late  fall  or  early  spring.  The  position  of  the  U.S.  portion  of 
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Table  6.  Load  distribution  ratio  on  ICS  at  selected  orientations. 


Position 

Test 

U.S. 

(N/m) 

ICS 

(Ibf/ft) 

Can.  ICS 
(N/m)  (Ibf/ft) 

Load  ratio 
U.S./Can. 

Orientation 

n 

1 

8 

645 

44.2 

794 

54.4 

0.81 

65 

1 

9 

727 

49.8 

910 

62.3 

0.80 

65 

2 

12 

1197 

82.0 

1543 

105.7 

0.78 

65 

3 

5a 

1466 

100.4 

1637 

112.1 

0.90 

85 

3 

Sb 

1515 

103.8 

1717 

117.6 

0.88 

85 

3 

6a 

1025 

70.2 

1143 

78.3 

0.90 

85 

3 

6b 

1048 

71.8 

1198 

82.1 

0.87 

85 

4 

7 

1231 

84.3 

1279 

87.6 

0.96 

105 

the  ICS  offered  some  flexibility  in  its  orientation  angle  with  respect  to  the  shoreline  and  its  uplake 
position. 

To  evaluate  the  U.S.  section  of  the  ice  control  structure,  the  Canadian  ICS  remained  stationary 
for  all  tests.  Since  it  is  very  difficult,  if  not  impossible,  to  reproduce  exactly  the  same  ice  floe 
positions  for  each  test,  the  shear  stress  on  the  ice  underside  varies  from  test  to  test  because  of  the 
different  floe  arrangements.  The  assumption  is  made  that  the  ice  conditions  are  generally  similar 
over  the  entire  model  and  that  any  change  in  conditions  for  each  test  is  also  similar  over  the  entire 
model.  This  assumption  allows  us  to  normalize  the  load  measured  at  the  U.S.  ICS  by  the  Canadian 
ICS  load.  Figure  19  shows  ttu  f.our  experimental  locations  for  the  ICS  orientation  test  and  Table 
6  gives  the  tabulated  results.  The  orientation  angle  is  measured  clockwise  from  the  shoreline  to  the 
ICS  as  was  shown  in  Figure  19. 

Figure  22  shows  quite  clearly  the  effect  of  the  orientation  angle  on  the  load  distribution  of  the 
U.S.  ICS.  The  smaller  the  included  angle  between  the  shoreline  and  the  ICS  (as  shown  by  positions 
1  and  2  in  Fig.  19),  the  lower  the  load  to  the  U.S.  ICS.  This  interpretation  implies  that  more  of 
the  load  is  being  taken  by  the  shoreline  along  the  U.S.  side.  A  practical  limit  is  reached  as  to  the 
minimum  included  angle,  because  the  ICS  must  increase  in  length  to  reach  the  same  stationary 
location  near  the  navigation  channel. 
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Figure  22.  Results  of  orientation  force  tests— ratio  of 
U.S.  I  Can.  loads  on  ICS. 
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The  conclusions  reached  from  these  tests  show  that  the  orientation  of  the  structure  should  have 
an  included  angle  of  slightly  less  than  90°  measured  clockwise  from  the  upstream  shoreline  to  the 
ICS  (i.e.  it  should  be  somewhere  between  position  1  and  position  3  in  Fig.  19). 

While  the  same  hydraulic  conditions  and  same  quantity  of  ice  were  maintained  and  the  only 
difference  between  tests  was  the  ice  floe  interaction,  the  resulting  drag  force  varied  by  a  factor  of 
roughly  1.5  for  four  tests  at  position  3.  Throughout  all  the  orientation  tests,  the  Canadian  ICS 
with  a  constant  position  took  loads  that  differed  by  a  factor  of  at  least  2.  These  few  tests  confirmed 
that  many  runs  had  to  be  undertaken  to  get  a  reasonable  range  for  the  force  distribution  in  mod¬ 
eling  as  well  as  in  the  field  data,  since  the  effective  shear  stress  appeared  to  be  highly  variable. 

Once  the  orientation  tests  were  completed,  the  location  of  the  ICS  was  set  according  to  the  gen¬ 
eral  alignment  of  position  3  in  Figure  19.  Table  7  represents  the  coordinates  of  the  122-m  (400- 
ft)  gap  using  the  Michigan  State  Survey  Coordinate  System.  The  other  gap  positions  may  be  cal¬ 
culated  readily,  as  the  gap  kept  the  same  bearing  and  its  center  coincided  with  the  centerline  of 
the  shipping  channel.  The  relative  bar  positions  1-1 1  are  given  in  Figure  23  for  all  tests  conducted. 
The  intermediate  locations  of  the  force  measuring  rods  for  tests  with  the  122  m  (400  ft)  gap  are 
also  given  in  Table  7. 

From  Table  7  the  positions  of  the  measuring  bars  for  the  other  tests  may  be  calculated.  The 
movement  of  these  rods  for  all  tests  would  be  less  than  ±  61  m  (200  ft)  from  their  location  as 
given  in  the  previous  table. 

The  difference  in  the  total  load  due  to  the  bars  in  slightly  different  locations  would  be  negli¬ 
gible.  Any  adjustment  in  the  load  distribution  would  be  small  due  to  the  small  change  in  velocity 
between  stream  lines  within  this61-m  (200-ft)  distance. 

Table  7.  Survey  coordinates  for  the  122-m  (400-ft)  gap. 


Location 

Michigan  state  coordinates 

Bar  position 

U.S.  (CS  Gap 

834  1  97 

553062 

4 

Can.  ICS  Gap 

834597 

553082 

5 

U.S.  Shoreline 

823751 

552575 

1 

Can.  Shoreline 

836175 

551086 

11 

833240 

552715 

2 

833738 

552922 

3 

848940 

552697 

6 

835244 

552204 

7 

835649 

551799 

8 

835639 

551761 

9 

836025 

551370 

10 

ANALYSIS  OF  ICE  DISCHARGE  DUE  TO  SHIP  TRANSITS 

The  most  critical  time  for  ice  discharge  through  an  opening  in  an  ice  control  structure  due  to 
ship  transits  is  when  the  ice  field  is  in  an  unconsolidated  state  (when  fragmented  randomly  shaped 
ice  floes  are  in  a  random  mixture).  The  planar  dimension  of  these  randomly  shaped  ice  floes  can 
vary  from  a  few  meters  to  hundreds  of  meters. 

The  analysis  of  the  ice  discharge  due  to  ship  transits  is  similar  in  some  respects  to  the  problem 
of  ice  arching  studied  by  Calkins  and  Ashton  (1975).  The  bleeding  of  ice  through  the  gap  in  an 
ice  boom  after  a  ship  transit  can  be  a  function  of  many  independent  variables: 

A r  -  f{o,  b,  k,  Vvbv  «M)  (13) 
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where  Ar  =  area  of  ice  released  per  vessel  transit 

a  =  average  floe  size  dimensions  from  eq  1 1 
b  =  gap  opening 

k  =  ice  type  and  ice  field  conditions 
V w  =  water  velocity 
<t>  =  vessel  direction 
l/y  =  vessel  speed 
bs  =  beam  of  ship 

{■  =  orientation  and  geometry  of  the  ICS. 

The  orientation  tests  performed  earlier  fixed  the  geometric  parameters  4/  and  %  at  constant 
values  for  all  the  further  experimental  work.  To  study  ice  discharge  it  was  decided  to  use  a  single 
vessel  of  fixed  beam  whose  dimensions  would  be  nearly  those  of  the  largest  ships  now  operating 
in  the  Great  Lakes.  The  water  velocity  (related  to  the  shear  stress  on  the  ice  cover)  was  relatively 
constant  for  the  flow  range  considered;  consequently,  all  tests  were  conducted  at  the  maximum 
flow  rate  and  maximum  water  levels  expected  at  Port  Huron. 

The  variables  can  be  reduced  to  a  nondimensional  grouping  without  bs,  Kw,  0,  p.  The  ice  dis¬ 
charge  normalized  by  the  gap  opening  yields 

Ar/b2  =  {alb,  VJVW,  k,  <I>,  b/bs)  .  (14) 

Three  different  fragmented  ice  types  k  were  evaluated:  freshwater,  urea-doped  and  polyethy¬ 
lene  plastic.  The  ship  direction  4>  was  either  upbound  or  downbound.  The  ratio  VJV^  varied 
over  a  small  range,  but  its  major  influence  was  how  the  ice  field  was  rearranged  by  the  ship  transits. 
The  last  parameter  in  eq  12  (b/bs)  was  found  to  be  insignificant  by  itself  in  analyzing  the  ice 
bleeding  data.  This  reduces  eq  4  and  the  ice  discharge  to  a  function  of  three  nondimensional 
quantities: 

Arlb2  =  f2(alb,k,<i>) .  (15) 

The  bleeding  of  ice  through  openings  in  ice  booms  due  to  ship  passages  has  been  documented 
for  the  St.  Marys  River  by  the  Detroit  District.  The  three  years  of  record  showed  an  average  sur¬ 
face  ice  discharge  of  0.35  b2  per  vessel  passage.  In  an  idealized  flume  study  of  ice  discharge  Calkins 
and  Ashton  (1975)  report  a  value  of  0.9  b 2  for  synthetic  square  ice  floes  without  regard  to  floe 
size  or  gap  opening.  Other  physical  model  studies  conducted  by  Acres  American  Inc.  (1975)  and 
Arctec  Inc.  (1978)  report  average  ice  discharge  values  of  6  b 2  and  2  b2,  respectively,  for  booms 
with  a  gap  of  69  m  (225  ft). 

The  effect  of  propeller  wash  on  ice  bleeding  was  not  investigated  because  its  effect  would  be 
minimal  due  to  the  limited  area  of  the  propeller  in  comparison  to  the  size  of  the  ice  floes.  How¬ 
ever,  another  concern  of  some  members  of  our  staff  was  that  the  modification  of  ship  stern  and 
the  alternating  direction  of  ship  passages  could  be  biasing  the  ice  bleeding  data.  To  evaluate  this 
possibility,  the  following  tests  were  conducted  with  vessel  transits  in  one  direction  (upstream) 
for  1)  the  ship  with  identical  bow  and  stern  sections  (tests  76  and  77),  and  2)  the  ship  with  normal 
bow  and  stern  sections  (test  78). 

The  tests  were  conducted  using  the  randomly  shaped  synthetic  ice  floes  with  an  ICS  gap  corre¬ 
sponding  91.5  m  (300  ft)  and  similar  hydraulic  conditions.  The  results  (Table  8)  show  that  there 
is  little  difference  in  ice  bleeding  for  upstream  transits  between  a  normal  ship  (test  78),  and  a  ship 
with  an  identical  bow  for  the  stern  (tests  76  and  77).  Tests  70,  74,  and  75  are  alternating  vessel 
passages  and  their  mean  ice  discharge  does  not  significantly  differ  from  that  of  tests  76,  77  and  78. 
Test  78  reflects  the  ship  with  a  normal  configuration  and  operated  in  one  vessel  direction  only.  The 
mean  discharge  value  lies  within  the  range  of  the  data  taken  for  tests  70,  74  and  75. 
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Table  8.  Test  parameters  for  plastic  ice  floes. 


Test 

Ice  thickness 
(m)  (ft) 

Floe  she  (a) 
(m)  (ft) 

Gap  opening  (b) 
(m)  (ft) 

a/b 

Remarks 

15 

0.54 

1.8 

34.5 

113 

189.2 

600 

0.188 

Square  floes 

16 

0.54 

1.8 

34.5 

113 

122 

400 

0.283 

17 

0.54 

1.8 

34.5 

113 

91.5 

300 

0.377 

Ice  rafting  noticed 

18 

0.54 

1.8 

34.5 

113 

198.2 

650 

0.174 

70 

0.54 

1.8 

17.4 

57 

91.5 

300 

0.190 

Randomly  stuped  floes 

71 

0.54 

1.8 

17.4 

57 

122 

400 

0.143 

72 

0.54 

1.8 

17.4 

57 

152.4 

500 

0.114 

73 

0.54 

1.8 

17.4 

57 

182.9 

600 

0.095 

74 

0.54 

1.8 

17.4 

57 

91 J 

300 

0.190 

75 

0.54 

1.8 

17.4 

57 

91.5 

300 

0.190 

76 

0.54 

1.8 

17.4 

57 

91.5 

300 

0.190 

All  upbound  transits 

77 

0.54 

1.8 

17.4 

57 

91.5 

300 

0.190 

All  upbound  transits 

78 

0.54 

1.8 

17.4 

57 

91.5 

300 

0.190 

Upbound  using  normal 

ship 


Table  9.  Test  parameters  for  laboratory  grown  fragmented 
ice  floes. 


Ice  thickness  Floe  size  (a)  Gap  opening  (b) 


Test 

(m) 

(W 

(m) 

(ft) 

(m) 

(ft) 

a/b 

24 

3.1 

19.7 

52 

122 

400 

0.129* 

25 

2.5 

26.8 

88 

122 

0.220f 

26 

1.02 

3.3 

27.7 

91 

122 

0.225 1 

27 

1.19 

3.9 

11.0 

36 

122 

400 

0.090f 

28 

0.72 

2.4 

16.2 

53 

91.5 

300 

0.1 76f 

29 

0.51 

1.7 

8.8 

29 

91.5 

300 

0.097 1 

30 

3.3 

12.2 

40 

91.5 

300 

0.132f 

31 

0.64 

2.1 

8.2 

27 

91.5 

300 

0.09  It 

32 

2.1 

8.5 

28 

91.5 

300 

0.093f 

33 

1.44 

4.7 

11.3 

37 

mm 

400 

0.092* 

34 

1.11 

3.6 

14.9 

49 

400 

0.121* 

35 

1.28 

4.2 

13.7 

45 

wSm 

400 

0.1 13f 

36 

0.51 

1.7 

15.5 

51 

400 

0.127 

37 

0.60 

2.0 

12.5 

41 

122 

400 

0.101 

39 

0.85 

2.8 

9.2 

30 

400 

40 

0.77 

2.5 

13.7 

45 

182.9 

600 

0.075 

41 

1.28 

4,2 

19.8 

65 

182.9 

600 

0.108 

42 

1.19 

3.9 

33.8 

111 

182.9 

600 

0.185 

60 

1.19 

3.9 

9.0 

29 

91.5 

300 

0.098** 

61 

0.68 

2.2 

39 

91.5 

300 

0.132** 

*  No  ice  rafting  during  tests, 
t  Ice  rafted  adjacent  to  ship  tuck. 
**  Urea  doped  ice 


All  the  data  for  the  ice  bleeding  tests  are  summarized  in  Tables  8  through  1 1  for  the  natural 
and  synthetic  ice  floes.  The  combined  tests  represent  over  1 000  vessel  transits  through  the  four 
ICS  openings.  In  the  tables,  the  ice  discharge  per  vessel  transit  (Ar)  has  been  normalized  by  the 
gap  opening  squared  (b2).  For  each  test  the  mean  and  standard  deviation  of  the  ice  discharge  per 
vessel  transit  in  one  direction  are  given.  Some  vessel  transits  produced  no  ice  discharge  and  the 
adjusted  values  of  ice  discharge  with  the  "zero  data”  removed  are  given  in  Tables  10  and  11.  Also 
computed  are  the  exceedance  levels  for  three  frequencies  of  90,  50  and  10%  from  a  fit  of  the 
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Table  10.  Ice  discharge  results  for  laboratory -grown  fragmented  ice  floes. 


Ice  discharge 

TP 


Adjusted 

Ar/b2 


Mean/std.  dev. 

90 % 

5096 

10% 

No.  of 
tests 

Skew 

coef. 

Mean/std.  di 

24U  0.49/0.45 

0.06 

0.35 

1.09 

11 

1.13 

0.49/0.45 

D  0.49/0.5 1 

0.05 

0.34 

1.13 

11 

1.43 

0.49/0.51 

25  U  0.14/0.19 

0.00 

0.06 

0.39 

19 

1.27 

0.19/0.20 

D  0.29/0.32 

0.00 

0.19 

0.70 

19 

1.23 

0.35/0.33 

26  U  0.04/0.06 

0.00 

0.03 

0.12 

37 

1.83 

0.07/0.06 

D  0.21/0.27 

0.00 

0.13 

0.51 

37 

2.44 

0.22/0.28 

27U  0.32/0.37 

0.02 

0.22 

0.74 

26 

1.78 

0.34/0.37 

D  0.54/0.71 

0.06 

0.36 

1.20 

26 

2.67 

0.54/0.71 

28U  0.06/0.J  1 

0.00 

0.03 

0.18 

34 

2.49 

0.90/0.12 

D  0.23/0.21 

0.04 

0.18 

0.48 

36 

1.47 

0.23/0.21 

29IJ  0.08/0.09 

0.00 

0.05 

0.22 

15 

0.18 

0.11/0.08 

D  0.13/0.07 

0.05 

0.11 

0.22 

IS 

0.39 

0.13/0.07 

30U  0.17/0.22 

0.01 

0.11 

0.40 

30 

2.36 

0.18/0.22 

D  0.37/0.39 

0.07 

0.29 

0.78 

30 

2.71 

0.37/0.39 

31 U  0.06/0.06 

0.00 

0.06 

0.14 

27 

0.55 

0.08/0.05 

D  0.21/0.12 

0.07 

0.18 

0.37 

27 

0.68 

0.21/0.12 

32U  0.16/0.21 

0.00 

0.11 

0.36 

30 

3.10 

0.18/0.21 

D  0.31/0.16 

0.14 

0.29 

0.51 

30 

1.49 

0.31/0.16 

33U  0.95/0.81 

0.18 

0.74 

1.94 

16 

1.26 

0.95/0.81 

D  1.31/0.65 

0.52 

1.18 

2.24 

16 

0.51 

1.31/0.65 

34U  0.55/0.36 

0.14 

0.46 

1.07 

20 

0.48 

0.55/0.36 

D  1.38/0.95 

0.52 

1.22 

2.42 

20 

2.36 

1.38/0.95 

35  U  0.23/0.34 

0.00 

0.12 

0.59 

29 

2.44 

0.25/0.35 

D  0.35/0.60 

0.01 

0.18 

0.88 

29 

2.80 

0.35/0.60 

36U  0.05/0.07 

0.00 

0.02 

0.13 

28 

1.26 

0.07/0.07 

D  0.14/0.23 

0.00 

0.08 

0.36 

28 

2.23 

0.15/0.23 

37U  0.20/0.23 

0.01 

0.13 

0.49 

22 

1.28 

0.21/0.23 

D  0.41/0.35 

0.04 

0.28 

0.97 

22 

0.29 

0.41/0.35 

39U  0.82/0.49 

0.31 

0.73 

1.43 

10 

0.87 

0.82/0.49 

D  0.80/0.57 

0.18 

0.65 

1.61 

10 

0.42 

0.80/0.57 

40U  0.32/0.2 3 

0.05 

0.25 

0.68 

10 

0.28 

0.32/0.23 

D  0.53/0.39 

0.12 

0.43 

1.04 

10 

0.87 

0.53/0.39 

41 U  0.34/0.18 

0.11 

0.29 

0.62 

5 

0.70 

0.34/0.18 

D  1.23/0.20 

- 

- 

- 

5 

- 

1.23/0.20 

42U  0.02/0.04 

0.00 

0.01 

0.06 

20 

2.21 

0.03/0.04 

D  0.36/0.36 

0.04 

0.25 

0.80 

20 

1.63 

0.36/0.36 

601)  0.26/0.50 

0.01 

0.14 

0.67 

15 

3.53 

0.26/0.50 

D  0.35/0.38 

0.02 

0.22 

0.84 

15 

1.02 

0.35/0.38 

61U  0.17/0.38 

0.00 

0.08 

0.44 

10 

3.13 

0.17/0.38 

D  0.19/0.31 

0.01 

0.11 

0.46 

10 

2.46 

0.19/0.31 

Sample 

size 


*  Static  forces  measured 
t  No  data  taken. 

**  Dynamic  forces  measured. 


individual  data  by  a  2  parameter  gamma  function  .  If  other  exceedance  levels  are  desired,  the  ice 
discharge  values  can  be  computed  using  the  normal  Pearson  Type  III  distribution  with  the  data 
from  Tables  10  and  11.  See  Appendix  C  for  further  explanation  of  the  data-fitting  procedure. 


Table  11.  Ice  discharge  results  for  plastic  ice  floes. 


Ice  discharge  _ Adjusted 


Meanjstd.  dev. 

9096 

A,/*2 

5096 

1096 

No.  of 

tests 

Skew 

coef. 

Ar/b* 

Meanjstd.  dev. 

Sample 

size 

15U  0.57/0.67 

0.04 

0.40 

1.32 

20 

2.j6 

0.60/0.68 

19* 

D  1.04/0.44 

0.39 

0.93 

1.82 

17 

-0.S9 

1.04/0.44 

17* 

16U  0.31/0.30 

0.00 

0.27 

0.69 

20 

1.30 

0.42/0.28 

15* 

D  0.94/0.43 

0.40 

0.86 

1.58 

20 

0.07 

0.94/0.43 

20* 

17U  0:19/0.30 

0.00 

0.00 

0.61 

21 

0.40 

0.44/0.31 

9* 

D  0.80/0.45 

0.24 

0.72 

1.46 

20 

0.26 

0.84/0.42 

19* 

18U  1.74/2.29 

0.08 

1.02 

4.36 

17 

1.50 

1.74/2.29 

17+ 

D  1.05/0.43 

0.59 

1.00 

1.56 

14 

1.30 

1.05/0.43 

14+ 

70U  0.07/0.12 

0.00 

0.00 

0.23 

25 

1.11 

0.16/0.15 

lit 

D  0.19/0.35 

0.00 

0.00 

0.61 

25 

1.39 

0.39/0.42 

12f 

71 SJ  0.15/0.17 

0.00 

0.09 

0.38 

25 

1.00 

0.19/0.18 

20f 

D  0.71/0.59 

0.11 

0.56 

1.48 

25 

1.58 

0.74/0.59 

24f 

72  U  0.64/0.49 

0.07 

0.45 

1.46 

20 

0.37 

0.64/0.49 

20f 

D  0.97/0.94 

0.22 

0.78 

1.92 

20 

2.53 

0.97/0.94 

20f 

73U  1.42/1.30 

0.10 

0.90 

3.40 

7 

0.68 

1.42/1.30 

7t 

D  1.71/1.13 

0.44 

1.42 

3.34 

7 

0.34 

1.71/1.13 

7t 

74U  0.04/0.08 

0.00 

0.00 

0.15 

10 

-0.59 

0.13/0.09 

3** 

D  0.25/0.33 

0.00 

0.15 

0.50 

10 

1.55 

0.27/0.34 

9** 

75  U  0.12/0.14 

0.00 

0.08 

0.30 

10 

0.47 

0.17/0.14 

7** 

D  0.81/0.38 

0.42 

0.76 

1.24 

10 

1.54 

0.81/0.38 

10** 

76U  0.13/0.11 

0.00 

0.11 

0.27 

20 

0.97 

0.15/0.10 

17f 

77U  <0.0 

- 

- 

- 

10 

- 

0.08 

It 

78U  0.10/0.10 

- 

0.075 

0.24 

20 

0.80 

0.14/0.09 

14+ 

*  Static  forces  measured, 
t  No  data  taken. 

**  Dynamic  forces  measured. 


Natural  ice 

The  ice  discharge  data  were  broken  down  into  upstream  and  downstream  vessel  transits.  Fig¬ 
ures  24  and  25  show  the  ice  discharge  (50%  exceedance  level)  as  a  function  of  ajb  for  downbound 
and  upbound  passages.  The  data  appear  to  be  scattered  at  first  glance  but  this  scatter  can  be  attrib¬ 
uted  to  ice  floes  of  larger  than  normal  dimensions.  Alongside  tests  where  the  ice  was  very  thick 
(>  1  m)  or  where  ice  floes  were  larger  than  20  m,  a  small  superscript  (f  or  1 )  is  indicated,  and  the 
superscript  data  points  appear  to  form  an  upper  bound  on  the  curves.  The  thicker  and  larger  floes 
did  not  raft  as  easily  as  the  thinner  and  smaller  sized  ice  floes.  This  rafting  of  the  ice  floes  enhanced 
the  re-arching  of  the  ice  floes  above  the  opening  after  each  transit,  thereby  minimizing  the  ice 
bleeding.  These  rafted  small  ice  floes  acted  as  a  very  cohesive  unit,  merely  separating  and  rejoining 
as  the  ship  passed  through  the  rafted  ice  field.  The  small  rafted  floes  effectively  acted  as  large  floes 
instead  of  individual  smaller  ones.  Figures  26a-c  show  a  sequential  view  of  a  downstream  vessel 
transit,  while  Figure  26d  is  a  closeup  of  the  rafted  ice  conditions  near  the  ship  track. 

The  effect  of  the  ice  field  conditions  on  ice  discharge  begins  to  disappear  when  the  alb  ratio  is 
less  than  0.1  of  the  upbound  transits.  The  data  are  more  greatly  scattered  for  the  downstream  tran¬ 
sits  but  the  same  trends  is  visible.  The  downstream  vessel  transits  caused  a  slightly  greater  ice  dis¬ 
charge  for  the  same  alb  ratio  than  the  upstream  transits  because  of  the  bow  push  of  ice  through 
the  opening  for  the  downbound  transits.  The  same  general  trends  are  also  apparent  using  the  mean 
ice  discharge  or  the  other  ice  discharge  exceedance  levels  (see  Fig.  27  for  upbound  ice  discharge  at 
the  10%  level).  Re-arching  of  the  ice  field  when  the  floe  size  to  gap  opening  ratio  alb  is  smaller 
than  0.075  appears  to  be  difficult  and  0.05  is  probably  the  lower  limit  for  ice  arching  to  be  success¬ 
ful  with  a  full  ice  field  behind  the  structure. 


35 


m  « 


mm**- 


-  *.«;K®** 


<-'■  ;..^%N,'-^,; 


»"■  ■  ■ 


9fel.,». 


$. .  «£' 


Figure  26.  Down  bound  ship  transit-fragmented  ice.  o-c)  Three  sequential  views  of  the  ship  transiting  downstream 
through  a  fragmented  ice  field  and  the  ICS.  Note  the  relatively  sharp-edged  ship  track,  d)  Close-up  of  the  rafted 
ice  floes  adjacent  to  the  ship  track. 


A  general  statement  regarding  the  ice  discharge  is  that  the  ice  discharge  occurring  at  50%  of  the 
time  for  an  unconsolidated  cover  will  not  exceed  ft2,  for  either  up  or  downbound  vessel  transits  if 
the  floe  size/gap  opening  ratio  is  greater  than  0.1.  If  the  ice  conditions  are  better  documented  with 
respect  to  flow  size,  condition  of  cover,  etc.,  then  a  refinement  in  the  ice  discharge  per  vessel  pas¬ 
sage  can  be  achieved  by  using  Figures  24  and  25. 


a/b 

Figure  2  7.  Upbound  ice  discharge—  1 0%  exceedance 
level— natural  Ice. 

Synthetic  ice 

Experiments  were  also  conducted  with  synthetic  ice  (polyethylene  plastic  with  a  specific  gravity 
of  0.92)  to  determine  if  it  could  serve  as  a  substitute  for  natural  ice.  Both  randomly  shaped  and 
square  uniform  floes  were  cut  from  the  large  sheets  of  polyethylene.  The  square  uniform  floes  of 
a  single  layer  would  be  re-arch  at  a/b  values  less  than  0.15  to  0.17;  however,  the  randomly  shaped 
floes  re-arched  at  alb  values  of  less  than  0.1.  Both  floes  had  the  same  prototype  ice  thickness  at 
0.54  m  (1.8  ft). 

The  square  uniform  floes  exhibited  a  much  higher  ice  discharge  for  the  same  a/b  ratio  than  the 
randomly  shaped  floes.  Figures  28  and  29  reflect  the  50%  exceedance  level  of  ice  discharge  for 
downstream  and  upstream  transits,  respectively,  while  Figure  30  uses  the  mean  ice  discharge  of 
upstream  transits.  The  mean  and  the  50%  exceedance  ice  discharges  reproduced  identical  graphs, 
displaced  only  by  the  skewness  in  the  data  at  each  alb  ratio.  The  ice  discharge  for  the  square  uni¬ 
form  floes  are  nearly  an  order  of  magnitude  greater  than  the  randomly  shaped  floes. 

The  same  trend  as  for  the  natural  random  floes  is  observed  for  the  random  synthetic  ice  floes: 
convergence  of  the  ice  discharge  rate  to  a  common  value  for  upstream  and  downstream  transits  as 
the  ratio  alb  approaches  0.05.  This  is  partially  due  to  the  large  gap  opening  that  allows  more  ice  to 
pass  through  the  gap  opening. 

Measurements  on  the  effective  depth  of  the  ice  rafting  near  the  ship  track  ranged  from  two  to 
four  floe  thicknesses  and  extended  in  scaled  width  up  to  60-90  m  (200-300  ft)  on  each  side  of  the 
ship  track  centerline.  This  rafted  zone,  if  greater  than  the  opening  in  the  ice  control  structure, 
definitely  retards  the  ice  bleeding  of  an  unconsolidated  ice  field  behind  the  structure.  If  the  open¬ 
ing  is  greater  than  120  to  180  m,  the  ice  rafting  due  to  ship  passages  is  not  as  important  a  factor 
for  reducing  the  ice  discharge. 

Figures  31a-b  show  the  ship  in  the  random  synthetic  ice  field  and  a  closeup  view  of  the  rafted 
ice  floes  adjacent  to  the  ship  track.  Figures  32a-d  show  sequential  photography  of  an  upbound 
vessel  transit  through  an  ice  field  composed  of  square,  uniform,  synthetic  ice  floes. 

The  ice  discharge  data  from  the  random  synthetic  ice  can  be  overlapped  with  the  data  from  the 
natural  ice  floe  test.  Figures  33  and  34  show  the  ice  discharge  per  vessel  transit  at  the  50%  exceedance 
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Figure  28.  Downbound  ice  discharge— 50%  exceedance 
level— synthetic  ice. 


(o)  Tests  70-751  Random 
(a)  Test*  76a77J  Shops* 

(•)  Tests  15-18!  Square 

and  Uniform 


Figure  30.  Upbound  ice  discharge— mean— synthetic  ice. 


a.  b. 

Figure  31.  Two  sequential  views  of  the  ship  traveling  upstream  through  the  ICS  and  the  rafting  of  the  plastic  ice. 
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Figure  33.  Downbound  ice  discharge,  50%  exceedance  level,  com¬ 
bined  natural  and  synthetic  random  floes  per  vessel  passage. 


Figure  34.  Upbound  ice  discharge,  50%  exceedance  level,  com¬ 
bined  natural  and  synthetic  random  floes  per  vessel  passage. 


level  for  downbound  and  upbound  vessel  transits  respectively.  The  random  synthetic  ice  data  are 
represented  by  the  solid  line,  which  splits  the  envelope  of  natural  ice  data  for  both  up  and  down- 
bound  transits.  This  similarity  is  very  significant.  In  model  studies  of  the  future  where  ice  discharge 
is  the  important  parameter  in  the  design  or  planning  of  a  hydraulic  structure,  synthetic  random¬ 
shaped  ice  floes  mav  be  used  with  assurance  that  they  represent  median  conditions  with  respect  to 
natural  random-shaped  ice  floes. 


FORCES  ON  THE  ICE  CONTROL  STRUCTURE 

The  forces  on  the  ICS  were  measured  in  two  ways,  depending  on  the  type  of  data  desired.  Static 
force  measurements  were  taken  at  the  1 1  bar  positions  prior  to  each  ship  transit  through  the  ice 
field  and  ICS  when  the  ice  cover  was  stationary  and  no  ice  was  flowing  through  the  opening.  The 
dynamic  measurements  of  forces  to  the  ICS  began  when  the  ship  approached  the  ice  field  and  con¬ 
tinued  for  a  period  of  either  33  or  50  seconds,  scanning  at  a  rate  of  3  samples/second  for  bars  1-9. 

Static  measurements 

The  static  measurements  were  taken  to  measure  the  changes  in  the  overall  force  levels  occurring 
after  each  ship  passage  through  the  ice  field  and  ICS.  Two  types  of  ice  were  used,  plastic  and  nat¬ 
ural  ice  floes,  and  the  data  were  evaluated  in  either  an  overall  force  analysis  on  the  structure  or  on 
an  individual  point  to  point  measurement  basis.  The  overall  force  method  resolves  the  forces  on 
each  of  the  two  ICS  structures  to  common  axes.  The  total  load  is  then  reduced  by  the  added  force 
from  the  weight  component  in  the  flow  direction  and  finally  divided  by  the  area  of  the  ice  cover. 

This  yields  an  effective  shear  stress  because  a  portion  of  the  downstream  ice  load  is  being  taken  by 
the  shoreline.  A  summary  of  die  data  for  some  of  the  tests  is  given  below  to  illustrate  the  low 
variability  within  a  test  without  physically  disturbing  the  ice  cover  in  the  total  redistributed  load 
to  the  ICS  for  varying  numbers  of  ship  passages  (Table  1 2). 

The  mean  effective  shear  stress  changes  from  one  test  to  another  because  re-arrangement  of  ice 
floes  creates  a  rougher  ice  underside  either  throughout  the  entire  model  or  quite  possibly  in  a  lo¬ 
calized  area  around  the  ICS  and  ship  track.  The  Student’s  t-test  done  on  any  two  pairs  of  data  may 
or  may  not  show  that  the  sample  means  are  from  the  same  population.  For  example,  runs  3  to  8, 

9  to  12  and  13  to  18  in  test  6  only  differ  in  that  the  ice  floes  were  repositioned  in  the  model  between 

Table  12.  Redistributed  total  force  on  ICS  (un¬ 
sealed  model  data) 


Test  numbers 


16-3-16-7 

16-3-16-8 

16-9-16-12 

16-13-16-18 

16- 19-16-25 

17- 2-17-7 
17-8-17-1 1 
24-5-24-8 

24- 10-24-15 

25- 3-25 -14 

25- 15-25-18 

26- 4-25-14 

27- 7-27-24 


Effective  shear 
stress  t  std,  dev. 
(N  ’ 


0.01 16  t  0.0008 

0.0110  ±  0.0015 
0.0268  ±  0.0016 
0.0086  t  0.0019 
0.0121  ±  0.0006 
0.0113  *  0.0001 
0.0139  ±  0,0022 
0.0126  ±  0.0020 
0.0133  ±  0.0019 
0.0074  t  0.0016 
0.0069  t  0.0014 
0.0120  ±  0.0022 
0.0265  ±  0.0021 


Vessel 

passages 


Ice  type 


Natural 

Natural 

Natural 

Natural 

Natural 

Natural 


tests  8  and  9,  and  1 2  and  1 3,  No  other  change  occurred,  yet  the  mean  effective  shear  stress  values 
differed  by  a  factor  of  2  to  3.  The  standard  deviations  were  always  less  than  20%  of  the  mean 
values  for  both  natural  ice  and  plastic  ice,  with  consistent  values  as  seen  in  the  tables.  If  the  shear 
stress  values  were  scaled  to  prototype  conditions,  the  range  for  the  means  would  be  0.587  to  2.28 
N/m2,  with  a  mean  of  1.14  N/m2,  only  very  slightly  lower  than  the  mean  value  measured  in  the 
field.  In  other  words,  die  measured  shear  stress  values  fall  within  die  range  of  the  calculated  field 
measurements. 

The  distribution  of  load  to  the  ICS  can  be  broken  down  in  terms  of  the  forces  taken  at  each  bar. 
The  following  tables  will  give  the  average  normal  load  at  the  measuring  points  in  terms  of  percentage 
of  the  total  normal  load.  The  total  normal  load  is  designated  as  the  sum  of  the  loads  normal  to  the 
orientation  to  each  of  the  two  ICS’s.  This  normal  load  component  closely  corresponds  to  the  direc¬ 
tion  of  flow  in  the  model,  as  the  ICS  was  purposely  situated  to  be  nearly  normal  to  the  flow  lines. 

In  fact  the  loads  could  also  be  approximated  in  a  radial  configuration  due  to  the  large  arc  that 
could  be  constructed  from  the  converging  flow  lines. 

Table  13  represents  the  average  percentage  of  normal  force  taken  at  each  measuring  position  for 
test  series  16, 17,  20,  22,  24,  25,  26,  and  27.  Positions  1  and  1 1  are  nearly  on  shore  and  carried  a 
very  small  percentage  of  the  load  because  the  ice  grounded  along  the  shoreline  and  transmitted  the 
load  to  the  shore  and  bed.  Position  bars  2  and  10  are  located  farther  into  the  lake  and  carried  loads 
equal  to  10%  and  6%  of  the  total,  respectively. 

The  two  rods  that  were  adjacent  to  the  ship  track  (4  and  5)  measured  roughly  10%  each  of  the 
total  load.  The  largest  loads  were  taken  by  rods  3  and  6  (21  and  16%)  which  were  one  position  away 
from  the  ends  of  the  U.S.  and  Canadian  ICS,  respectively.  This  was  not  unexpected  as  rods  3  and 
6  have  more  effective  ice  boom  lengths  than  rods  4  and  5.  The  remaining  two  rods  along  the 
Canadian  ICS  measured  13%  each  of  the  total  load.  Overall  42%  of  the  mean  load  was  taken  by 
the  U.S.  ICS  and  58%  by  the  Canadian  ICS,  but  variations  would  occur  within  each  test.  The  max¬ 
imum  percentage  measured  for  one  side  of  the  ICS  within  one  of  the  runs  was  58%  taken  by  the 
U.S.  side,  and  70%  for  1  run  by  the  Canadian  ICS.  No  distinction  was  made  between  upbound  and 
downbound  vessel  transits  for  these  data. 

These  percentages  can  be  further  analyzed  at  a  local  level  by  considering  the  maximum  percent¬ 
age  taken  by  any  one  measuring  position  during  all  of  the  1 12  vessel  transits  for  these  8  tests.  Table 
14  gives  the  maximum  percentage  taken  by  any  rod  during  all  of  these  runs.  Position  3  during  one 
of  the  10  runs  in  test  22  took  55%  of  the  total  measured  load  for  that  run.  The  next  highest  per¬ 
centage  was  34%  taken  by  position  8-9.  The  remaining  positions,  excluding  the  two  end  rods,  had 
maximum  percentage  loads  ranging  from  22%  to  32%.  These  loads  represent  redistributed  static 
loads  on  the  structure  after  a  vessel  has  transited  through  the  ICS. 


Table  13.  Average  percentage  of  total  normal  force  taken  at  each 
position. 

No,  of 


Test 

7 

2 

3 

4 

5 

6 

7 

8-9 

10 

77 

runs 

16 

3 

14 

16 

9 

8 

8 

12 

12 

17 

0 

17 

5 

10 

18 

17 

4 

13 

8 

18 

7 

0 

23 

20 

0 

5 

22 

9 

9 

15 

13 

22 

6 

0 

8 

22 

2 

9 

23 

14 

5 

10 

18 

13 

7 

0 

24 

0 

7 

19 

14 

S 

20 

10 

3 

0 

11 

25 

2 

13 

30 

13 

6 

18 

12 

3 

2 

0 

16 

26 

1 

10 

25 

19 

6 

16 

21 

2 

0 

11 

27 

D 

KEj 

12 

8 

14 

24 

12 

11 

8 

0 

23 

10  21  11 
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Table  14.  Maximum  percentage  of  normal  load  taken  after 
any  run  within  each  test.  (The  maximum  value  in  each 
column  is  italicized.) 


Test 

1 

2 

3 

4 

5 

6 

7 

8-9 

10 

11 

16 

8 

16 

22 

13 

13 

12 

19 

22 

22 

0 

17 

11 

16 

21 

19 

9 

19 

10 

21 

12 

0 

20 

0 

7 

35 

15 

13 

20 

15 

24 

9 

0 

22 

5 

15 

55 

23 

10 

18 

31 

26 

13 

0 

24 

0 

10 

23 

21 

12 

30 

18 

34 

4 

0 

25 

3 

26 

49 

27 

10 

25 

21 

7 

10 

0 

26 

2 

15 

38 

24 

17 

28 

32 

5 

0 

0 

27 

1 

15 

17 

16 

25 

32 

15 

14 

10 

0 

Dynamic  force  measurements 

As  mentioned  earlier,  continuous  data  were  taken  on  the  forces  generated  by  a  ship  moving 
through  the  ice  field  and  the  ICS  opening.  Tests  were  conducted  with  fragmented  floes  of  plastic, 
freshwater  ice  and  doped  urea  ice.  Only  the  urea-doped  ice  tests  will  be  described,  as  the  flexual 
strength  was  scaled  according  to  values  of  typical  freshwater  ice,  6  to  7  x  106  Pa.  Comparison 
between  the  ice  types  has  not  been  completed. 

Table  15  is  a  summary  of  the  dynamic  data  taken  during  test  61  for  20  vessel  transits,  10  up- 
bound  and  10  downbound.  Each  run  consists  of  measuring  the  two-directional  force  150  times 
for  each  position  (1-9)  during  a  50-s  time  frame.  The  vessel  could  traverse  up  or  down  through 
the  ice  field  and  ICS  in  about  20  s  depending  on  the  ship  speed  setting.  The  values  in  Table  14 
represent  the  ratio  of  maximum  load  recorded  to  the  mean  load  for  the  10  runs  in  both  directions; 
i.e.  the  sample  population  is  1 500  readings  per  bar  per  direction.  The  x  direction  represents  the 
normal  load  registered  to  the  ICS  while  the  y  direction  reflects  the  traverse  loading  at  each  point. 
Due  to  the  relatively  symmetrical  pattern  offered  by  the  ICS  and  its  relation  to  the  flowlines,  the 
y-direction  loads  were  considerably  lower  than  the  x-direction  loads. 

The  ratios  are  generally  consistent  except  for  the  y  direction  at  position  3  and  both  x  and  y 
directions  at  position  9.  This  basic  data  show  that  the  primary  reason  for  the  large  ratios  at  posi¬ 
tions  3  and  9  are  the  relatively  low  mean  values.  Measured  values  for  some  of  the  runs  constituted 
loads  of  less  than  2  g  (or  less  than  2700  Ibf  prototype). 


Table  15.  Ratio  of  peak  load  to  mean  load  for  urea- 
doped  ice  at  each  measuring  position.  Data  were 
taken  during  test  61  for  10  upbound  and  10  down- 
bound  vessel  transits. 

1  0  upbound  1 0  downbound 


X 

Y 

X 

Y 

Max/mean 

Max/mean 

Max/mean 

Max/mean 

Bar 

ratio 

ratio 

ratio 

ratio 

1 

1.4 

2.6 

1.8 

4.9 

2 

4.2 

3.8 

3.8 

4.2 

3 

2  2 

>10 

1.6 

>10 

4 

1.9 

3.6 

3.4 

4.7 

5 

2.3 

2.6 

1.9 

1.8 

6 

1.9 

2.2 

1.8 

2.2 

7 

1.7 

1.9 

1.5 

1.8 

8 

2.4 

2.9 

2.1 

2.0 

9 

>10 

9.3 

>10 

>10 

45 


The  two  major  conclusions  drawn  from  Table  14  are  1)  there  is  little  difference  in  the  ratio  of 
the  peak  to  mean  loads  for  upbound  or  downbound  ship  passages  at  each  position  and  2)  the  ratios 
also  appear  to  be  relatively  consistent  from  position  to  position. 

These  data  suggest  that  the  vessel-ice-floating  structure  interaction  is  independent  of  vessel  direc¬ 
tion.  The  rearrangement  of  the  ice  floes  caused  by  the  vessel  is  the  dominant  dynamic  loading 
feature  with  a  fragmented  ice  cover.  This  same  conclusion  (loads  applied  to  the  ice  boom  are  inde¬ 
pendent  of  vessel  direction)  was  observed  for  the  St  Marys  River  ICS  by  Perham  (1978).  The  dy¬ 
namic  load  to  the  ICS  in  terms  of  its  static  load  offered  by  the  ice  cover  should  average  on  the 
order  to  3  to  5  times  the  static  load  conditions,  but  extreme  values  of  greater  than  10  were  noted. 

The  information  given  in  Tables  12-14  may  be  applied  in  determining  the  load  on  the  ICS  due 
to  the  ice  cover  water  drag  and  ship  interaction  with  the  fragmented  ice  cover.  Once  the  loading 
to  the  structure  on  a  force  per  linear  dimension  basis  (i.e.  Ibf/ft  or  N/m)  has  been  determined,  the 
load  distribution  can  be  computed  from  Tables  12  and  13.  When  the  load  distribution  has  been 
calculated,  the  dynamic  load  factors  from  Table  14  can  be  applied  to  each  position  to  account  for 
ship-ice-structure  interaction  at  each  loading  position.  A  factor  of  5  would  appear  appropriate  for 
all  positions. 

The  effect  of  a  strong  wind  on  the  force  distribution  of  the  ICS  with  a  ship  running  through 
a  solid  ice  cover  was  not  considered  in  the  physical  model,  although  inferences  on  these  effects  can 
be  provided.  The  ship-ice  floe  interaction  is  a  very  localized  phenomenon  and  any  additional  wind 
stress  would  be  acting  over  a  very  small  area  of  floes  disturbed  by  the  ship  passage.  It  is  our  belief 
that  the  effect  of  wind  on  the  interaction  between  the  ship  and  the  ice  cover  with  respect  to  the 
dynamic  loading  would  be  minimal  at  the  ICS  due  to  the  very  local  conditions. 


POTENTIAL  ADDITIONAL  SHEAR  STRESSES 

Since  the  size  of  area  to  be  modeled  was  restricted  by  the  building  dimensions,  the  influence 
of  the  other  portions  of  the  lake  on  the  shear  stress  was  evaluated  mathematically.  A  logarithmic 
velocity  profile  distribution  was  used  with  the  maximum  velocity  taken  at  0.75 Du  where  Dx  is  the 
depth  of  flow.  The  roughness  height  kt  was  taken  as  1.5  cm  (0.05  ft).  The  velocity  profile  data 
taken  at  Port  Huron  suggested  these  values  (which  are  conservative). 

An  idealized  situation  of  flow  in  the  lake  was  mathematically  represented.  The  discharge  was 
assumed  constant  at  5667  m3/s  (200,000  ft3/s),  and  the  flow  depth  increased  linearly  from  6  m 
(20  ft)  at  1.6  km  (1  mile)  into  the  lake  to  a  depth  of  15  m  (50  ft)  at  a  distance  of  29  km  (18  mile) 
into  the  lake.  The  effective  width  of  the  proposed  ICS  from  shoreline  to  shoreline  was  estimated 
at  1 220  m  (4000  ft).  The  shear  stress  was  calculated  by  using  eq  5  over  an  incremental  lake  length 
of  AL  =  1 140  m  with  an  effective  width  of  1 20  m  (see  Fig.  35).  The  downstream  force  F\,  the 
product  of  the  shear  stress  on  the  ice  cover  r,  and  the  ice  cover  area  A  j,  is  summed  for  each  suc¬ 
cessive  ice  cover  length  (/  increment)  extending  into  the  lake: 

n 

F,=  2  r!  A[  .  (14) 

/-I 

The  calculation  was  halted  when 


<  0.01  . 


There  are  no  field  data  available  at  this  time  to  confirm  the  assumptions  made  on  the  velocity 
or  shear  stresses  this  far  into  the  lake.  The  purpose  for  this  calculation  was  to  estimate  the  relative 
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LAKE  HURON 


Figure  35.  Schematic  view  of  additional  water  drag  compo~ 
nents. 


(Ibf/f)  (N/m) 


let  Cover 


Figure  36.  Estimated  load  to  ICS  with  surface  area  not  included 
in  the  mode I  limits. 


magnitude  of  die  possible  additional  load  due  to  an  ice  cover  by  the  water  shear  stresses  outside  the 
limits  of  the  model  area. 

Figure  36  shows  the  results  of  using  eq  14  with  the  assumptions  given  earlier.  The  additional 
load  decreases  quite  rapidly  with  distance  away  from  the  ICS;  the  1%  level  is  calculated  to  be  ap¬ 
proximately  5.69  km  (18.675  ft)  outside  the  limit  of  the  model.  This  additional  load  is  relatively 
small,  only  48.8  N/m  (36  Ibf/ft)  for  an  ICS  with  a  length  of  1220  m  (4,000  ft).  This  calculation 
assumes  that  no  load  is  taken  by  an  arch  or  through  shear  along  the  ice/ice  interface  boundary. 

This  additional  load  is  an  upper  bound  on  the  additional  load  due  to  water  stresses  not  considered 
in  the  model  due  to  space  limitations. 
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ANTICIPATED  ICE  CONDITIONS  WITH  ICS 


The  ICS  will  undoubtedly  initiate  an  ice  cover  as  soon  as  the  first  ice  appears  at  Port  Huron. 

The  ICS  will  change  die  local  ice  conditions  along  the  U.S.  and  Canadian  shorelines  so  that  1 )  the 
ice  will  form  earlier  and  the  thickness  of  shore  ice  pileup  may  increase,  2)  the  distribution  and 
thickness  of  shore  ice  pileup  will  be  altered,  and  3)  the  extension  of  the  present  shear  walls  and 
shore-fast  ice  near  the  inlet  to  the  river  should  move  into  the  lake  and  then  taper  back  to  their 
present  location  within  1  to  2  km  upstream  of  the  ICS  along  each  shoreline  (see  Fig.  8).  The  un¬ 
rafted  zone  of  ice  beyond  the  shorefast  ice  would  be  free  to  move  uplake  if  a  sufficient  shear  stress 
were  available.  The  shorefast  ice  would  not  be  moved  uplake  due  to  its  anchored  resistance  along 
the  shoreline  and  bed.  The  returning  ice  from  uplake  would  then  strike  and  arch  between  the 
shorefast  ice.  The  shoreline,  being  the  most  rigid  component  of  the  shorefast  ice,  ICS  and  shore¬ 
line  system,  would  absorb  the  majority  of  the  load. 

To  minimize  forces  on  the  ICS,  the  shorefast  ice  should  remain  in  place  to  serve  as  the  buffer 
strip  between  large  ice  floes  arriving  from  uplake  and  the  ICS.  The  present  observation  of  ice  along 
the  shoreline  confirms  the  hypothesis  that  the  shorefast  ice  remains  until  the  very  end  of  the  season. 
With  an  ICS  restricting  its  movement  downstream,  the  ice  will  remain  in  place  until  it  melts  or  is 
destroyed  by  open  water  wave  action. 

The  ship  passages  through  the  ICS  will  affect  the  ice  conditions  upstream.  The  ship  track  will 
become  clogged  with  brash  ice  at  times.  Open  leads  will  also  develop  behind  the  ICS,  both  within 
the  ship  track  and  in  other  areas.  These  open  leads,  possibly  extending  up  to  1  km  in  length,  will 
be  confined  by  the  shear  walls  of  the  shorefast  ice.  The  open  leads  will  form  by  repeated  ship 
traffic,  causing  small  amounts  of  ice  to  pass  through  the  opening,  or  possibly  from  thermal  melting 
of  the  ice  underside. 

Downstream  ice  conditions  in  the  St.  Clair  River  will  be  altered  by  prohibiting  ice  discharge 
from  Lake  Huron.  Assuming  an  average  thickness  of  15  cm  (6  in.)  over  the  ice  season  and  the 
average  seasonal  ice  discharge  of  900  km2,  about  1.35  x  108  m3  of  ice  would  be  held  back  by 
the  ICS. 


CONCLUSIONS 

The  results  of  this  modeling  effort  have  led  to  the  following  general  and  specific  conclusions  re¬ 
garding  modeling  of  an  ICS  in  a  subzero  environment. 

1.  Special  construction  considerations  have  to  be  used  for  the  bed  of  the  model.  Use  of  fiber¬ 
glass  resin  only  or  swimming  pool  paint  over  the  conventional  sand-base  mortar  skin  will  work  well. 

2.  In  the  open  water  calibration,  downstream  adjustments  at  the  gate  did  not  affect  upstream 
flow  conditions  near  the  ICS  or  navigation  channel  and  only  at  the  inlet  were  adjustments  possible. 
Flow  acceleration  toward  the  opening  was  dominant,  and  the  flowlines  were  found  to  be  very  sim¬ 
ilar  for  all  flows  evaluated. 

3.  The  ice  cover  calibration  was  achieved  by  obtaining  a  properly  scaled  shear  stress  which  was 
close  to  the  mean  of  data  collected  in  the  field.  The  many  tests  showed  a  wide  variation  in  the 
shear  stress  values,  which  were  calculated  from  direct  measurements  of  force  and  divided  by  the 
ice  surface  area.  Both  plastic  and  natural  ice  offered  the  same  effective  shear  stress,  both  having 
the  same  absolute  range  of  values,  from  roughly  0.007  N/m2  to  0.027  N/m2.  Repetitive  testing 
was  therefore  found  to  be  an  absolute  necessity. 

4.  For  the  random-shaped  natural  floes,  the  ice  conditions  upstream  of  the  structure,  the  floe 
size/gap  opening  ratio  (alb),  and  the  vessel  direction  were  the  predominant  independent  variables 
affecting  the  ice  discharge  per  vessel  transit  (Ar/b2). 
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a.  The  surface  conditions  of  the  ice  floes  upstream  of  the  ICS  (i.e.  whether  the  floes  were  rafted 
or  not)  was  very  dominant,  as  a  non-rafted  ice  field  exhibited  two-  to  five-fold  greater  ice 
discharge  per  vessel  transit  than  a  rafted  ice  field  for  similar  hydraulic  and  ICS  conditions. 

b.  The  ratio  of  floe  size  to  gap  opening a/b  was  very  significant  in  determining  the  ice  discharge. 

A  change  in  this  ratio  from  0.1  to  0.2  showed  a  10-fold  decrease  in  the  ice  discharge  per  vessel 
transit,  regardless  of  direction.  When  the  ratio  alb  was  greater  than  0.25,  the  discharge  was 
very  minimal  (<0.10  for  the  50%  exceedance  level).  The  minimum  a/b  ratio  for  re-establishing 
a  stable  arch  with  rafted  floes  was  approximately  0.075.  The  non-rafted  floes  approached  a 
minimum  ratio  value  of  0.10. 

c.  For  the  non-rafted  ice  fields,  downstream  vessel  transits  had  generally  higher  ice  discharges 
than  upstream  transits,  primarily  due  to  the  bow  push  of  ice  through  the  opening. 

d.  The  thicker  and  larger  ice  floes  used  in  these  tests  did  not  raft  as  easily  as  the  thinner  anu 
smaller  sized  floes.  These  tests  indicated  floes  with  scaled  sizes  greater  than  1  m  thick  and 

20  m  to  a  side  did  not  raft  as  quickly,  if  at  all.  When  the  width  of  the  rafted  ice  zone  adjacent 
to  each  side  of  the  ship  track  was  less  than  the  width  of  the  gap  opening,  the  effects  of  the 
rafted  floes  in  reducing  the  amount  of  ice  discharged  were  not  always  significant  The  scaled 
width  of  rafted  ice  was  generally  60  to  90  m  (200  to  300  ft)  from  each  side  of  the  ship  track 
centerline,  for  a  120-  to  180-m  total  effective  width. 

5.  The  two  synthetic  ice  floe  types  (random  and  square  shaped)  displayed  large  differences  in 
ice  discharge  for  similar  flow  and  gap  opening  conditions  with  respect  to  the  vessel  direction  and 
the  flow  size  to  gap  opening  ratio. 

a.  The  downstream  vessel  transits  exhibited  higher  ice  discharge  per  vessel  passages  than  upstream 
passages  for  both  square  and  randomly  shaped  floes.  The  downstream  transits  through  the 
randomly  shaped  floes  showed  a  50%  increase  in  ice  discharge  over  the  upstream  passages. 

The  randomly  shaped  synthetic  floes  exhibited  ice  discharge  values  similar  to  those  of  the 
randomly  shaped  natural  ice  floes.  Both  upstream  and  downstream  curves  for  the  synthetic 
random  floes  lie  within  the  bed  of  data  for  the  natural  randomly  shaped  ice  floes. 

b.  The  square  uniform  size  ice  floes  exhibited  a  significantly  higher  ice  discharge  by  a  factor  of 
10  than  the  randomly  shaped  floes.  The  downstream  ice  discharge  for  square  floes  was  en¬ 
tirely  controlled  by  bow  push  when  alb  >  0.20. 

c.  Square  uniform  synthetic  floes  are  significantly  different  in  their  ice  discharge  characteristics, 
not  rafting  as  easily  and  not  being  able  to  re-arch  at  a/b  <  0.1 5.  Random  synthetic  floes 
would  re-arch  at  alb  ~  0.1,  a  value  slightly  higher  than  that  for  natural  floes. 

6.  The  forces  to  the  ice  control  structure  from  the  ice  sheet  varied  considerably  from  test  to 
test.  Both  plastic  ice  floes  and  natural  ice  floes  exhibited  the  same  range  of  effective  shear  stress 
values. 

a.  The  average  percentage  of  the  total  normal  load  taken  by  each  of  the  measuring  bars,  exclu¬ 
sive  of  the  two  end  points,  was  generally  1 0  to  1 5%,  except  for  one  bar  (3)  which  took  roughly 
20%.  This  bar  was  situated  in  the  model  where  the  water  velocity  was  very  high  compared 

to  the  other  bars. 

b.  The  maximum  normal  load,  as  a  percentage  of  the  total  load  taken  by  any  one  bar  during 
the  112  tests,  was  55%  by  bar  3.  The  other  remaining  eight  bars  (2, 4,  5, 6,  7,  8,  9, 10)  took 
between  22  and  34%  of  the  total  load  during  any  one  test  run. 

c.  The  dynamic  loads  registered  by  the  1 1  bars  in  the  x-direction  (i.e.  parallel  to  the  flowlines) 
was  three  to  five  times  the  average  load  for  a  particular  run.  The  level  of  the  dynamic  loads 
measured  was  independent  of  the  vessel  direction.  Downstream  vessel  transits  did  not  create 
higher  dynamic  loads. 

7.  Future  model  studies  that  need  to  know  only  the  ice  discharge  rates  through  surface  openings 
may  use  synthetic  ice  cut  into  randomly  shaped  sizes.  The  polyethylene  plastic  of  0.92  specific 
density  is  a  good  substitute  for  natural  ice  for  ice  discharge  studies. 
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8.  The  model  results  for  average  ice  discharge  per  vessel  transit  are  consistent  with  the  field  data  J 

obtained  from  the  St  Marys  River.  The  only  field  data  on  ice  discharge  averaged  0.35  b 2  per  vessel  \ 

passage  in  the  St.  Marys  River  while  the  physical  model  ranged  from  0.1  to  1.0  b2  depending  on  | 

the  ratio  alb,  ice  type,  ice  conditions  upstream  of  the  boom  and  vessel  direction.  This  close  agree¬ 
ment  with  the  field  data  gives  confidence  that  die  modeling  effort  was  successful  for  both  randomly  j 

shaped  natural  and  synthetic  ice  floes. 

9.  The  orientation  of  any  structure  placed  at  Port  Huron  should  be  similar  to  the  testing  orien-  ] 

tation,  normal  to  the  flowlines.  The  gap  opening  of  1 22  m  (400  ft)  works  well  and  it  offers  a  larger 

margin  of  safety  (from  ship  damage)  than  a  61-m-  (20-ft-)  wide  opening.  1 
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APPENDIX  A.  APPLICATION  OF  MODEL  RESULTS 


This  section  is  written  for  the  user  of  the  results  derived  from  the  two  model  studies  and  the 
other  non-modeling  efforts.  The  majority  of  the  data  was  purposely  non-dimensional ized  for  easier 
interpretation  at  this  stage  in  the  report  The  detailed  design  of  the  ICS  will  be  done  by  another 
agency,  but  the  basic  design  criteria,  loads  and  ice  discharge  through  the  ICS  are  contained  in  this 
section. 

The  opening  in  the  ICS  will  allow  for  vessels  and  ice  to  pass  through  the  gap.  If  the  ice  field  con¬ 
ditions,  vessel  direction  and  gap  opening  are  known,  the  amount  of  ice  passing  during  vessel  transits 
can  be  calculated.  The  ice  cover  conditions  at  Port  Huron  can  be  highly  variable,  from  a  solid  cover 
one  day  to  fragmented  broken  floes  the  next.  The  worst  condition  for  ice  bleeding  under  natural 
and  ship-disturbed  conditions  would  be  the  prevalence  of  small  fragmented  ice  floes  whose  size  to 
gap  openir.g  ratio  is  less  than  0.1 .  These  particular  conditions  are  difficult  to  control,  but  fortun¬ 
ately  they  occur  only  a  small  percentage  of  the  time.  Currently,  a  solid  ice  cover  of  Port  Huron 
does  not  always  release  ice  downriver  when  disturbed  by  ship  transits.  The  photographic  documen¬ 
tation  brings  this  point  to  light. 

If  the  ice  conditions  at  Port  Huron  can  be  forecast,  observed  or  measured  similarly  to  the  docu¬ 
mentation  of  ice  conditions  in  the  model,  detailed  use  of  the  model  ice  discharge  figures  is  justified. 
If  the  ice  conditions  are  too  variable,  then  a  conservative  approach  can  be  used  to  estimate  ice  dis¬ 
charge  through  the  opening. 

Example  1 :  A  decision  was  reached  to  design  for  a  gap  opening  of  1 22  m  (400  ft).  Determine 
the  amount  of  ice  bleeding  for  100  round-trip  vessel  passages. 

Procedure 

1.  Determine  the  ice  field  conditions  for  the  winter  season  in  terms  of  floe  size. 

2.  Determine  the  number  of  ship  passages  for  each  floe  size  category. 

3.  Determine  the  50%  level  of  ice  discharge  per  vessel  transit  and  multiply  by  the  number  of 
upbound  and  downbound  transits  for  each  floe  size  to  gap  opening  ratio. 

Solution  1 

1.  Assume  three  ice  field  conditions 

a.  20  days,  a  =  15  m— non-rafted  January 

b.  40  days,  a  >  30  m— rafted  January-February 

c.  30  days,  a  =  20  m— non-rafted  March 

2.  One  hundred  round-trip  vessel  transits 

a.  40  January 

b.  20  January-February 

c.  40  March 

3.  Use  the  upper  and  lower  boundaries  in  Figures  24  and  25  for  non-rafted  and  rafted  ice  fields, 
respectively. 

a.  alb  =  0.1 3— non-rafted 

Ar!b 2  *  10  and  Ar/b 2  •<  0.5  (up  and  down) 

A-  40  x  (122)2  [1.0 +  0.5]  =0.89  km2 

b.  alb  =  0.25 

/4r/62  0.02  and  0.05  (up  and  down) 

A r  20  x  (122)2  [0.02  +  0.05]  =  0.02  km2 

c.  alb  =  0.16 

Arlb 2  0.5  and  0.2  (up  and  down) 

A,  40  x  (122)2  [o  s  +  o.2]  =  0.42  km2 

Total  1.33  km2 
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Solution  2 

One  could  use  a  maximum  value  of  approximately  1.5  for  the  ice  discharge  regardless  of  vessel 
direction  and  calculate  an  upper  bound  condition. 

Ice  released  =  100  x  (122)2  [1.5  +  1.5]  =  4.46  km2  (1.7  mile2) 

Example  2:  What  would  be  the  ice  discharge  through  the  gap  opening  under  natural  conditions 
without  ship  transits  and  before  the  ice  arched? 

A  solid  ice  cover  could  not  be  considered,  only  the  fragmented  cover.  The  ice  floe  size  to  gap 
opening  would  have  to  be  small,  i.e.  <  0.1,  and  a  very  low  ice  concentration  of  floes  at  the  gap. 

This  condition  would  be  most  prevalent  during  freezeup  or  after  an  ice  breakup  due  to  wind  action. 

A  conservative  time  duration  over  the  entire  winter  season  might  be  20  days,  and  10%  of  the  time 
the  ice  concentration  might  be  as  great  as  0.5  before  arching.  Using  a  transport  water  velocity  of 
0.6  m/s  (1.85  ft/s)  through  the  gap 

A  f  A  Cl  iAA  0.6  IT1  86400  S 

An  =  20  x  (10%)  [0.5]  122  m  x  — - —  x  —  — — 

An  =6.3  km2  (2.5  mile2). 

Example  3:  What  effects  would  wind  have  on  the  ice  bleeding  under  natural  conditions  and 
during  ship  transits? 

1.  Ice  bleeding  under  natural  conditions  will  continue  until  the  ice  discharge  or  ice  concentra¬ 
tion  at  the  gap  opening  for  a  given  a/b  ratio  exceeds  a  critical  value  as  determined  by  Calkins  and 
Ashton  (1975).  The  velocity  and  direction  of  the  moving  ice  floe  is  dependent  on  the  shear  stress 
offered  by  the  water  and/or  wind. 

The  ice  concentration  at  the  gap  is  independent  of  the  ice  floe  velocity  or  direction.  Consequently 
wind  will  have  little  direct  effect.  If  waves  in  art  open-water  environment  are  set  up  by  the  wind, 
the  ice  will  tend  to  move  along  the  opposite  shore  or  ice  edge.  This  would  impose  a  nonuniform  ice 
concentration  at  the  gap  opening,  and  the  re-arching  conditions  would  be  unfavorable.  Only  when 
the  ice  reacts  with  both  ice  edges  along  the  opening  will  the  arch  form. 

2.  Evaluating  the  effect  of  wind  on  the  ice  discharge  during  ship  transits  follows  the  same  line 
of  reasoning  given  above.  The  disturbed  ice  field  is  very  small  and  re-arching  of  the  ice  floes  is  de¬ 
pendent  upon  the  floe  size/gap  opening  ratio,  floe  conditions  and  to  some  extent  the  vessel  direction. 

Example  4:  What  are  the  design  loads  to  the  ICS  from  the  water-induced  forces  and  how  are  they 
distributed? 

1.  The  ICS  is  sited  normal  to  the  direction  of  the  flow  and  this  is  the  optimum  condition.  The 
shear  stress  developed  on  the  ice  sheet  resulted  in  loads  of  650  N/m  (45  ibf/ft)  to  1825  N/m  (125 
Ibf/ft)  of  ice  control  structure  length.  The  additional  load  to  the  structure  for  the  area  not  included 
in  the  model  could  reach  another  450  N/m  (30  Ibf/ft).  This  gives  a  total  static  load  of  roughly  2190 
N/m  (150  Ibf/ft)  of  ice  control  structure  length. 

2.  The  two  ice  control  structures  (U.S.  and  Canadian)  are  465  and  755  m  long  respectively  for 
a  total  length  of  1 220  m. 

a.  Z.T  =  total  load  =  1220  m  x  2190  N/m  =  2672  kN 
U.S.  =  43%  (£t)  =1149  kN 

Canadian  =  57%  (£.T)  =  1523  kN. 

b.  The  ratio  of  each  ICS  length  to  the  total  length  is  another  check  on  the  ICS  load  distribu¬ 
tion,  if  a  relatively  uniform  load  distribution  from  the  ice  fields  is  assumed. 

465 

U.S.  =  j22o  x  2671.8  =  1018  kN 
Canadian  =  -j~-  x  2671.8  =  1654  kN  . 

The  load  distribution  is  consistent  for  both  methods. 
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Example  5:  The  location  of  the  ICS  will  be  moved  to  a  different  position  than  that  tested  in 
the  model.  What  corrections  to  the  data  need  to  be  made? 

1.  The  ICS  is  in  a  minimum  downstream  location  along  the  U.S.  shoreline.  It  would  be  possible 
to  shift  the  Canadian  ICS  downstream  until  a  velocity  of  0.6  to  0.7  m/s  (2.0-2.3  ft/s)  is  reached. 

2.  The  change  in  the  loading  to  the  ICS  can  be  calculated  by  multiplying  the  average  shear  stress 
over  the  change  in  surface  area  to  account  for  inclusion  or  omission  of  extra  load.  The  increase  in 
load  assuming  a  200-m  position  movement  downstream  is  roughly  100  N/m,  i.e.  less  than  5%  of 
the  total  used  in  example  4. 

3.  The  ice  bleeding  would  be  barely  affected  unless  the  submergence  velocity  in  the  gap  region 
was  exceeded  for  a  particular  floe  size.  An  upstream  or  downstream  movement  of  the  ICS  will  not 
affect  the  reported  data  by  more  than  5-10%,  if  that 

Example  6:  How  will  Tables  13, 14  and  1 S  assist  in  the  design  of  an  ICS? 

1.  Table  1 3  can  be  used  regardless  of  the  type  of  ICS  selected.  It  gives  the  percentage  of  the 
total  normal  load  (provided  the  structure  is  situated  normal  to  the  fiowlines)  at  each  measuring 
point  Since  the  structure  was  physically  divided  at  the  gap  opening,  the  percentage  of  the  total 
load  taken  by  the  U.S.  and  Canadian  sections  can  also  be  found  in  Table  13.  The  U.S.  section  con¬ 
tained  bars  1-4,  while  the  Canadian  portion  had  the  remainder  (5-1 1 ).  A  range  in  the  average  data 
can  be  extracted  from  Table  1 3. 

2.  Table  14  reflects  the  maximum  percentage  of  the  normal  load  taken  at  any  one  time  at  a 
particular  measuring  bar  during  one  of  the  1 1 2  transits.  The  maximum  value  has  been  highlighted. 
These  data  reflect  some  of  the  localized  ice  interaction  phenomena  very  near  the  gap  opening  as  a 
result  of  the  ship  passages.  The  data  in  Table  14  can  be  used  with  a  floating  ice  control  structure 
concept  but  not  directly  if  a  rigid  structure  is  planned. 

3.  Tables  14  and  15  can  be  used  if  a  floating  ice  control  structure  is  considered.  The  position 
of  the  measuring  bars  can  be  assumed  to  be  the  end  points  of  the  individual  floating  sections,  as 
they  were  in  the  model.  An  individual  floating  section  might  be  as  long  as  120  to  190  m  composed 
of  smaller  floating  elements  such  as  timbers  or  pontoons  to  retain  the  ice. 

Both  Tables  14  and  1 5  are  measures  of  the  extreme  loads  taken  at  any  one  point  This  percent¬ 
age  can  be  compared  in  Table  13  to  the  average  value  measured  for  that  measuring  bar.  The  maxi¬ 
mum  value  is  generally  two  to  three  times  the  mean  value  for  all  bars. 

Table  15  is  slightly  different  than  Table  14.  The  data  contained  in  Table  15  are  the  ratio  of 
maximum  loads  to  the  mean  load  at  each  position  for  both  x  and  y  directions  using  urea-doped 
ice.  These  data  were  abstracted  from  the  dynamic  force  measurements  during  ship  transits  through 
the  ICS  and  ice  field.  The  ratios,  except  for  the  two  end  bars  and  a y  direction  on  bar  3,  are  con¬ 
sistent  with  die  values  one  would  calculate  by  considering  the  data  in  Table  13  and  14.  There  are 
exceptions  where  the  ratio  is  close  to  four,  but  generally  the  peak  load  to  average  load  ratio  is  2 
to  3. 

4.  The  design  of  a  rigid  structure  to  retain  the  ice  at  Port  Huron  would  be  based  on  the  crushing 
strength  of  the  ice;  this  load  of  1380  to  2760  kN/m  (200-400  lbf/in.2)  far  exceeds  the  loading  con¬ 
sidered  in  this  study  using  fragmented  ice  floes  contained  by  a  floating  structure. 
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APPENDIX  B.  SUGGESTED  ADDITIONAL  STUDIES 


The  following  section  discusses  areas  where  additional  data  collection  is  necessary  to  substantiate 
some  findings  in  the  physical  model  section.  Also,  areas  are  suggested  where  more  research  is  needed 
to  verify  some  of  the  assumptions  in  the  wind  stress  analysis. 

Physical  model 

The  calibration  of  the  physical  model  for  ice-covered  conditions  requires  knowledge  of  the  shear 
stress  exerted  by  the  water  to  the  underside  of  the  ice  sheet.  Since  only  10  measurements  were 
available,  additional  data  need  to  be  gathered  to  verify  the  shear  stress  values  obtained  in  February 
1979.  The  decay  in  the  shear  stress  as  one  goes  farther  into  the  lake  would  shed  light  on  the  assump¬ 
tions  used  in  the  analysis  of  the  influence  of  the  portions  of  the  lake  outside  the  study  area  that 
could  contribute  load  to  an  ice  control  structure. 

To  minimize  any  doubt  about  the  open  water  calibration,  another  drogue  study  for  a  lower  dis¬ 
charge  (i.e.  <  180,000  cfs)  should  be  undertaken.  More  drogue  lines,  including  additional  lines 
closer  to  the  two  shorelines,  should  be  used. 

Winter  conditions 

The  available  data  on  ice  movement  into  the  St  Clair  River  are  probably  the  most  comprehen¬ 
sive  ever  taken  to  date  at  one  site.  To  fully  understand  the  ice  processes  at  Port  Huron,  the  following 
observations  and  measurements  would  provide  the  necessary  backup  information  for  the  work  al¬ 
ready  completed  on  ice  arching  and  for  evaluating  potential  environmental  effects:  1)  limits  of  the 
ice  arch  location,  2)  shoreline  ice  pileup  documentation,  3)  shoreline  profiles  of  the  bed,  4)  meteoro¬ 
logical  data  near  the  site,  5)  in-situ  ice  stress  measurements,  and  6)  critical  pressure  determinations 
for  ice  deformation. 

1.  The  ice  arch  location  throughout  the  entire  ice  season  would  help  verify  the  wind  stress 
modeling  work.  Determining  the  ice  conditions  that  lead  to  stabilizing  a  series  of  ice  arches,  partic¬ 
ularly  if  ship  traffic  is  occurring,  would  require  an  almost  daily  observation  by  air  and  close  coordi¬ 
nation  of  ship  movements. 

2.  The  shoreline  ice  pileup  occurs  now  without  a  structure  in  place.  Thus,  existing  ice  pileups 
and  windrows  need  to  be  documented  in  case  severe  ice  conditions  arise  with  an  ice  boom  in  place. 
Physical  measurements  of  these  pileups  and  their  distribution  over  distance  and  time  are  needed. 

The  distribution  of  these  events  can  be  mapped  for  aerial  observations  or  photographic  documen¬ 
tation.  The  physical  measurements  can  be  conducted  by  ground  surveys  with  the  proper  ice  augers, 
but  a  remote-sensing  radar  device  may  offer  more  flexibility  and  greater  coverage.  CRREL  person¬ 
nel  have  conducted  one  set  of  shore  pileup  measurements  at  Port  Huron  and  Sarnia  with  the  im¬ 
pulse  radar.  It  appears  from  the  initial  playback  of  data  that  the  technique  of  data  collection  is 
very  critical,  as  height,  speed  and  possibly  the  grounded  ice  are  important  in  gathering  reliable  data. 
The  radar  method  looks  promising  but  data  collection  techniques  need  to  be  refined.  Processing 
the  data  and  interpreting  the  output  will  require  field  verification. 

3.  The  shoreline  bed  profiles  will  establish  the  baseline  data  on  the  processes  taking  place  over 
a  winter  season,  as  well  as  during  the  entire  season.  The  ice  pileup  conditions  probably  disrupt 
the  bed  conditions,  but  how  much  and  how  frequently  are  not  known.  To  answer  these  questions 
would  require  ground  surveys  prior  to  and  following  the  ice  season  at  several  cross  sections  with 
tight  vertical  and  horizontal  control  at  each  section.  One  midwinter  observation  might  be  necessary 
if  ice  conditions  are  very  erratic.  Both  shorelines  would  need  to  be  observed  and  a  good  set  of  low- 
level  vertical  photographs  taken  every  two  years  to  quantify  any  noticeable  shoreline  retreat  over 

a  5-  to  10-mile  segment. 

4.  The  measurement  of  in-situ  ice  stress  and  how  it  correlates  with  the  prevailing  wind  direction 
and  speed  is  a  link  to  verify  portions  of  the  modeling  concept  The  details  of  how  to  make  these 
measurements  and  where  to  make  them  should  be  worked  out 
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5.  Determination  of  the  critical  pressure  that  initiates  the  ice  cover  to  deform  is  needed.  As 
mentioned  in  the  wind  stress  modeling  section  this  critical  stress  needs  to  be  determined  for  ice 
ridging,  rafting  and  shore  pileup  conditions.  A  laboratory  study  in  a  large  refrigerated  basin  with 
an  instrumented  towing  carriage  and  side  walls  instrumented  to  determine  the  total  force  balance 
would  be  necessary. 


APPENDIX  C.  DERIVATION  OF  ICE  DISCHARGE  FREQUENCY  ANALYSIS 


Ice  discharge  in  square  kilometers  per  vessel  transit  through  the  ICS  opening  could  vary  from 
zero  to  some  upper  boundary  value.  Theoretically,  the  upper  boundary  value  could  be  unrestrained, 
i.e.  infinite.  From  the  data  range  one  could  automatically  eliminate  a  normal  distribution  function. 
Since  the  raw  data  exhibited  a  definite  right  skewness  some  form  of  an  exponential  function  is  sug¬ 
gested.  Hydrologists  frequently  use  the  two-parameter  gamma  frequency  distribution  function  for 
data  analysis  of  extreme  events:  flood  flows,  rainfall,  etc. 

Due  to  the  noticeable  right-skewness  in  the  data  and  the  universal  acceptance  of  the  gamma  dis¬ 
tribution  in  statistical  hydrology,  it  was  assumed  to  be  a  reasonable  distribution  for  fitting  the 
data.  The  two-parameter  gamma  function  is  also  called  the  Pearson  Type  ill,  and  tables  are  found 
in  most  hydrology  texts  dealing  with  the  use  of  tables  for  the  incomplete  gamma  function  based  on 
frequency  factors. 

The  2-parameter  gamma  functions  has  the  following  distribution: 


fly)  =  — —7  y 
0“r(«) 


t-l  *>->-/<?  &>  0 

a>  0 


where  y  =  random  variable  ice  discharge  value 
0  =  scaling  parameter  ofy 
a  =  shape  parameter 
T  =  usual  gamma  function. 

The  estimation  of  parameters  a  and  0  was  evaluated  by  two  methods.  The  maximum  likelihood 
estimator  (MLE)  technique  advanced  by  Thom  (1957)  was  found  to  be  superior  over  the  method 
of  moments  and  used  in  this  analysis.  The  MLE  for  the  gamma  function  is 

1  +4/3  (In  Y-  1  liny) 

Q= - _  |  - - 

4  (In  Y  -  ^  Sin  Y) 

where  0=  Y/a 

Y-  mean  value 
N=  number  of  observations 

The  method  of  moment  estimation  for  a  is  Y2/S 2,  where  5*  is  the  variance  of  they  variable  and 
0  is  calculated  in  the  same  manner. 

The  integral  of  the  gamma  probability  density  function  yields  the  incomplete  gamma  function, 


1  («» 0;  Y)  =  J 
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where  Y is  the  value  of  y  at  the  upper  limit  of  the  integration.  Defining  f  =  y/0;0rfr  =  dy  and 
changing  the  variables  to  functions  of  t  yields 

V70  ttm  . 

dt 

Jo 

with  further  simplification  to 


57 


m 


P(a,  m=-^  f 
d 

P(a,m  =  T^  y(a,m. 

The  chi-square  distribution  is  also  related  to  the  incomplete  gamma  function  where 
P  (a,  YIP)  *  P  (x2|  v)  “  7  (a,  V70) 

X2  =  2K//3 


v  =  2a 


and  since  P  (x2|  can  be  approximated  by  a  series  expansion  for  small  values  of  a 
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The  evaluation  of  the  series  was  performed  by  factoring  the  recurrence  interval  until  the  abso¬ 
lute  difference  between  two  successive  computations  was  less  than  0.0001.  Stability  problems  did 
occur  with  large  values  of  a  (>20)  and  when  the  area  under  the  density  function  approached  0.9999. 
Also,  a  uniform  data  set  was  encountered  where  a  ~  35;  the  series  could  not  handle  this  alpha  value. 

The  gamma  function  T(x)  for  0  <  x  <  1,  was  also  evaluated  from  a  series  expansion.  The  form 
was  a  polynomial  approximation 


T  (x+1 )  =  x\  =  1  +  a1  x  +  a2x2  +  <73x3  +  o4x4  +  a5xs 

with  «71  =  -0.5748646 
a2  =  0.9512363 
o3  =  -0.6998588 
o4  =  0.4245549 
as  =-0.1010678 


The  values  of  the  gamma  function  for  n  >  1  were  computed  using  the  recurrence  relation 

T  (n  +z)  =  (n  -1  +*)  [n  -  2  +z) ...  (1  +z)  T(1  +z) 

where  T  (1  +z)  is  computed  from  the  previous  polynomial  approximation  for  values  less  than  1. 
Hand  calculations  for  these  quantities  are  unthinkable;  consequently  a  computer  program  for  the 
Hewlett-Packard  9845B  was  written  to  evaluate  the  above  series  and  values  of  the  probabilities  at 
the  desired  levels. 

The  ultimate  goal  was  to  determine  the  values  of  ice  discharge  per  vessel  transit  at  the  various 
levels  of  exceedance.  The  ice  discharge  data  were  separated  into  two  categories  of  vessel  direction 
and  ice  type. 

The  computed  values  of  means,  standard  deviations,  skew  coefficients,  etc.,  can  be  found  in 
Tables  10  and  1 1.  The  method  of  moments  estimation  was  evaluated  and  the  discharge  values  were 
generally  close  to  those  computed  by  the  MLE  method.  All  exceedance  level  discharge  values  re¬ 
flect  the  MLE  technique. 


A  straightforward  method  to  use  is  the  "frequency  factor"  approach  when  the  data  are  already 
tabulated.  This  is  obtained  through  a  mathematical  exercise  of  the  term  Yy  =  py  and  one  arrives  at 

k7  =  k+sJa- 

where  Yy  =  values  of  ice  discharge  at  exceedance  level  y{Ty) 

Y-  mean  ice  discharge 
5^  =  variance  of  ice  discharge 
t  -a 

K  =  frequency  factor  =  — —  =  f  (Ty,  Cs) 

Cs  =  sample  skew  coefficient. 


Tables  of  K  values  for  the  Pearson  Type  III  function  can  be  found  in  most  recent  hydrology  texts. 
To  find  K,  the  skew  coefficient,  and  probability  levels  are  needed.  The  arithmetic  skew  coefficient 
for  all  data  is  found  in  the  tables. 

There  are  several  data  bases  that  contained  zero  ice  discharge  for  vessel  transits.  The  statistical 
quantities  were  computed  for  the  non-zero  events  in  the  data  base  and  then  adjusted  according  to 
the  relation 


Pa  («.  j)=pt  +*t  ( P(«,  mi 


where  ( Y/(})  =  adjusted  probability  for  a  value  <  Y/P 


s  1  =  1- 


N 


n0  =  number  of  event  with  zero  ice  discharge 
N  =  total  number  of  ice  discharge  data  points  including  zeros. 


The  mathematical  functions  for  the  incomplete  gamma  function  were  obtained  from  the  U.S. 
Dept  of  Commerce  { 1 970). 


A  facsimile  catalog  card  in  Library  of  Congress  MARC 
format  is  reproduced  below. 
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